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IV 


ABSTRACT 


The  IEEE  802.1  la  standard  uses  a  eoded  orthogonal  frequeney  division  multi¬ 
plexing  (COEDM)  scheme  in  the  5-GHz  band  to  support  data  rates  up  to  54  Mbps.  The 
COEDM  was  chosen  because  of  its  robustness  to  multipath  fading  affects.  In  the  stan¬ 
dard,  convolutional  codes  are  used  for  error  correction.  This  thesis  examines  the  perform¬ 
ance  of  the  COEDM  system  with  variable  rate  Reed-Muller  (RM)  error  correction  codes 
with  a  goal  to  reduce  the  peak-to-average  power  ratio  (PAPR).  Contrary  to  the  expecta¬ 
tions,  RM  codes  did  not  provide  expected  improvement  in  PAPR  reduction.  Peak  clip¬ 
ping  and  Hanning  windowing  techniques  were  investigated  in  order  to  reduce  the  PAPR. 
The  results  indicate  that  a  tradeoff  exists  between  the  PAPR  and  the  bit-error  rate  (BER) 
performance.  Although  peak  clipping  yielded  considerable  reduction  in  PAPR,  it  required 
high  signal-to-noise  ratios.  On  the  other  hand,  Hanning  windowing  provided  only  a  small 
reduction  in  PAPR  with  reasonable  BER  performance. 
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EXECUTIVE  SUMMARY 


Since,  wireless  local  area  networks  (WLAN)  were  increasingly  becoming  more 
and  more  important  for  broadband  wireless  communieations  in  both  military  and  com- 
mereial  application,  the  main  goal  of  this  thesis  was  to  investigate  the  performanee  of  a 
coded  orthogonal  frequency  multiplexing  (COFDM)  system  utilizing  Reed-Muller  (RM) 
codes.  In  order  to  aehieve  this  goal,  three  different  types  of  RM  codes  and  for  eomparison 
a  rate  one-half  eonvolutional  code  are  investigated  through  simulations.  Also  Hanning 
windowing  and  peak  clipping  techniques  were  investigated  for  reducing  the  peak-to- 
average  ratio  (PAPR)  in  COFDM  systems. 

The  performance  of  the  COFDM  system  was  studied  under  indoor  and  oudoor 
channel  environments  using  five  different  channels.  The  bit-error  performance  eurves  for 
eaeh  of  the  RM  codes  and  convolutional  code  were  obtained. 

The  results  showed  that  the  COFDM  system  is  robust  in  indoor  ehannel  environ¬ 
ments.  However,  for  the  outdoor  channel  environments,  the  system  requires  higher  sig- 
nal-to-noise  ratios  (SNR)  to  aehieve  the  same  bit-error  rate  (BER)  performance  as  in  the 
indoor  case.  Also,  BER  performance  curves  showed  that  the  RM  codes  provide  better 
performanee  in  indoor  channel  environment  at  low  SNR.  However,  convolutional  codes 
provide  better  performanee  in  outdoor  channel  environments. 

Reed-Muller  eodes  are  straightforward  to  implement,  and  they  provide  a  wide 
range  of  coding  options.  However,  RM  codes  did  not  provide  the  expected  improvement 
in  PAPR. 

The  addition  of  Hanning  windowing  and  peak  elipping  improve  PAPR  reduction. 
While  improvement  with  Hanning  windowing  is  limited,  peak  elipping  provides  remark¬ 
able  reduetion  in  PAPR  but  at  the  eost  of  inereased  bit  error. 
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I.  INTRODUCTION 


Wireless  loeal  area  networks  (WLAN)  are  inereasingly  beeoming  more  and  more 
important  for  broadband  wireless  eommunieations  in  both  military  and  eommereial  appli- 
eations.  The  IEEE  802.1  la  WLAN  standard  uses  the  orthogonal  frequeney  division  mul¬ 
tiplexing  (OEDM)  to  support  data  rates  from  6  to  54  Mbps  and  operates  in  the  5-GHz 
band.  OEDM  is  a  multiearrier  transmission  method  in  whieh  a  high-rate  data  sequenee  is 
split  into  multiple  low  rate  subbloeks,  and  eaeh  subbloek  of  data  is  modulated  onto  sepa¬ 
rate  subearriers.  Speeifieally,  IEEE  802.1  la  uses  eoded  OEDM  (COEDM)  beeause  of  its 
robustness  under  multipath  fading  eonditions. 


A,  OBJECTIVE 

It  is  well  known  that  multiearrier  eommunieation  systems  suffer  from  high  peak- 
to-average  power  ratio  (PAPR)  beeause  of  the  eonstruetive  interference  of  the  signals. 
There  are  numerous  methods  to  reduce  the  PAPR  of  the  COEDM  signal.  One  of  the  miti¬ 
gation  techniques  is  to  use  Reed-Muller  (RM)  codes. 

The  main  goal  of  this  thesis  was  to  analyze  the  COEDM  system  utilizing  RM 
codes.  In  order  to  achieve  this  goal,  three  different  types  of  RM  codes  and  for  comparison 
a  rate  one-half  convolutional  code  were  investigated  through  simulation  studies.  Also, 
Hanning  windowing  and  peak  clipping  techniques  were  investigated  for  reducing  the 
PAPR  in  OFDM  systems. 

Although  binary  pulse  shift  keying  (BPSK),  quadrature  pulse  shift  keying 
(QPSK),  16  quadrature  amplitude  modulation  (16-QAM),  and  64-QAM  techniques  are 
used  in  IEEE  802.1  la  WLAN  systems,  only  the  QPSK  modulation  technique  was  imple¬ 
mented  in  the  simulations  here.  Using  four  different  simulation  steps,  the  performance  of 
the  COEDM  system  was  studied  under  both  indoor  and  outdoor  channel  environments. 
The  bit-error  rate  (BER)  performance  curves  for  each  of  the  RM  codes  and  the  convolu¬ 
tional  code  were  obtained. 
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B,  RELATED  WORK 

The  RM  coding  algorithms  in  [1]  provided  a  starting  point  for  the  implementation 
of  encoding  and  decoding  algorithms  of  RM  codes  in  the  simulations  here.  The  Hanning 
windowing  and  the  peak  clipping  algoritms  used  in  this  work  were  first  introduced  in  [2], 

The  MATLAB  simulation  code  used  in  this  thesis  was  first  developed  by 
Roderick  [3]  for  line-of-sight  ship-to-ship,  ship-to-shore  and  ship-to-relay  type  of  con¬ 
nectivity.  Tan  [4]  modified  the  MATLAB  code  in  [3]  to  develop  a  simulation  of  the  IEEE 
802.1  la  standard’s  physical  layer  for  indoor  channel  environments.  Convolutional  codes 
with  hard  decision  decoding  technique  were  used.  The  author  used  the  MATLAB  code  in 
[4]  as  a  starting  point  for  developing  the  simulations  in  this  thesis. 


C.  THESIS  ORGANIZATION 

Chapter  II  presents  the  basics  of  COEDM  and  its  schematic  structure,  the  defini¬ 
tion  of  PAPR,  and  a  general  overview  of  the  various  functional  blocks  in  an  OEDM  sys¬ 
tem.  Interleaving  and  guard  interval  are  explained  as  part  of  the  COEDM  system.  Chapter 
III  focuses  on  the  coding  scheme.  The  Hanning  windowing  and  peak  clipping  techniques 
are  described.  The  simulation  methodology  and  test  results  are  presented  in  Chapter  IV. 
Einally,  the  thesis  concludes  with  Chapter  V. 

Appendix  A  describes  the  MATLAB  architecture.  The  transmitter,  receiver  and 
channel  parts  of  the  simulation  model  are  explained,  and  MATLAB  functions  are  de¬ 
scribed.  Appendix  B  includes  the  MATLAB  code. 
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II.  CODED  ORTHOGONAL  FREQUENCY  DIVISION 
MULTIPLEXING  (COFDM)  AND  PEAK-TO-AVERAGE  POWER 

RATIO  (PAPR) 


This  chapter  presents  the  basics  of  the  eoded  orthogonal  frequency  division  mul¬ 
tiplexing  (COFDM)  teehnique  and  the  advantages  and  disadvantages  of  COFDM  eom- 
munieation  systems.  As  it  is  one  of  the  main  disadvantages  of  the  multicarrier  communi- 
eation  systems,  the  peak-to-average  power  ratio  (PAPR)  in  OFDM  systems  is  also  dis¬ 
cussed. 

A,  CODED  ORTHOGONAL  FREQUENCY  DIVISION  MULTIPLEXING 

(COFDM) 

In  COFDM,  the  total  bandwidth  is  divided  into  N  nonoverlapping  frequeney  sub- 
ehannels.  Eaeh  subehannel  is  modulated  with  a  different  symbol  and  then  the  N  subehan- 
nels  are  frequeney-multiplexed.  While  eonventional  frequeney  division  multiplexing 
(FDM)  sueeessfully  avoids  speetral  overlap  of  channels  and  eliminates  interehannel  inter- 
ferenee  (ICI),  it  makes  the  frequeney  speetrum  usage  ineffieient.  OFDM  uses  the  spee- 
trum  more  efficiently  than  the  eonventional  frequeney  division  multiplexing  (FDM)  by 
allowing  frequency  overlapping. 

Figure  1  illustrates  the  differenee  between  the  nonoverlapping  and  overlapping 
multiearrier  modulation  teehniques.  It  ean  be  seen  that  almost  half  of  the  bandwidth  is 
saved  by  overlapping  the  speetra. 
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OFDM 


FDM 


Figure  1.  Frequency  Efficiency  of  OFDM  over  FDM  Illustrated  for  Two  and  Three 

Subchannels  (After  Ref.  5.) 

OFDM  also  requires  that  the  subchannel  carrier  frequencies  are  orthogonal.  Two 
subcarriers  are  orthogonal  if  their  inner  product  is  zero.  The  COFDM  subcarrier  signals 
can  be  arranged  such  that  the  sidebands  of  the  carriers  overlap  but  the  signals  can  still  be 
received  without  ICI.  In  COFDM,  cross  talk  among  the  subcarriers  is  prevented  by  mak¬ 
ing  the  subcarriers  orthogonal.  In  other  words,  each  subcarrier  must  be  spaced  at  intervals 
of  1  /  Tj  where  is  the  symbol  duration  for  each  subcarrier.  A  representative  COFDM 

signal  is  shown  in  Figure  2.  Individual  COFDM  subchannel  response  is  shown  in  Figure 
2(a)  and  the  COFDM  signal  spectrum  for  A  =  4  is  shown  in  Figure  2(b).  It  can  be  seen 
that  the  spectra  of  the  subcarriers  are  not  separated  but  overlapped  in  frequency. 

The  discrete  Fourier  transform  (DFT)  is  used  to  accomplish  frequency  multiplex¬ 
ing  and  subcarrier  orthogonality.  The  transmitter  uses  an  IDFT  operation  to  modulate 
subblocks  of  data  symbols.  A  DFT  is  used  at  the  receiver  where  the  signal  is  sampled  at 
the  center  frequency  of  each  subchannel,  and  the  transmitted  data  are  recovered  without 
ICI.  From  Figure  2,  it  can  be  seen  that,  at  the  center  frequency  of  a  subchannel,  the 
sidelobes  due  to  all  the  other  channels  produce  a  null. 
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Figure  2.  Spectra  of  an  OFDM  Subchannel  and  OFDM  Signal  (From  Ref  5.) 

The  fast  Fourier  transform  (FFT)  algorithm  is  used  to  implement  both  the  IDFT  at 
the  transmitter  and  the  DFT  in  the  receiver,  thus  making  it  computationally  efficient. 

The  COFDM  divides  a  symbol  sequence  into  N  symbol  groups  so  that  each  group 
contains  a  symbol  sequence  of  rate  1/N  of  the  original  sequence.  The  COFDM  signal 
consists  of  N  orthogonal  subcarrier  signals,  where  each  of  subcarrier  is  modulated  by  a 
different  symbol  sequence.  Lengthening  the  symbol  duration  helps  eliminate  the  use  of 
an  equalizer  in  the  receiver. 

COFDM  reduces  the  effects  of  frequency  selective  fading  or  narrowband  interfer¬ 
ence.  If  fading  or  interference  takes  place  in  a  single  carrier  system,  it  can  fail  the  whole 
system  while  only  some  of  the  subcarriers  will  be  affected  in  a  multicarrier  system.  Addi¬ 
tionally,  by  using  error  correction  coding,  the  COFDM  system  performance  can  be  fur¬ 
ther  improved. 

A  COFDM  system  uses  the  frequency  spectrum  very  efficiently  since  it  allows 
overlapping.  In  COFDM,  the  use  of  cyclic  extension  prevents  intersymbol  interference 
(ISI)  and  ICI  problems.  Dividing  the  channel  into  subchannels  makes  the  COFDM  sys¬ 
tem  more  robust  to  frequency  selective  fading  effects  while  channel  coding  and  interleav¬ 
ing  make  it  more  robust  to  errors  due  to  the  channel  noise.  These  features  of  the  COFDM 
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are  considered  the  key  advantages.  However,  besides  these  advantages,  the  system  has 
also  drawbacks  as  it  requires  power  amplifiers  with  a  high  PAPR  in  order  to  provide  a 
large  effective  range  and  is  prone  to  frequency  offset  relative  to  single  carrier  systems. 

The  efficiency  in  frequency  spectrum  usage  and  the  robustness  to  multipath  fad¬ 
ing  make  COFDM  desirable  in  a  number  of  applications,  such  as  digital  audio  broadcast¬ 
ing  (DAB),  high  rate  digital  subscriber  line  (HDSL),  very  high  rate  digital  subscriber  line 
(VHDSL),  asymmetric  DSL  (ADSL),  HDTV  terrestrial,  IEEE  802.11  and  HiperLAN/2, 
general  switched  telephone  network  (GSTN),  cellular  radio,  and  digital  video  broadcast¬ 
ing  (DVB-T).  [6] 


B,  COFDM  SCHEMATIC  DIAGRAM 

The  basic  principle  of  COEDM,  as  mentioned  in  the  previous  section,  is  to  divide 
a  high-rate  data  stream  into  N  lower  rate  streams  and  to  transmit  them  at  the  same  time 
over  a  number  of  subcarriers.  Since  the  symbol  duration  is  increased,  the  relative  amount 
of  dispersion  in  time  caused  by  multipath  delay  spread  is  decreased.  Intersymbol  interfer¬ 
ence  (ISI)  is  another  problem,  which  can  almost  be  eliminated  by  introducing  a  guard 
time  in  every  COEDM  symbol.  In  order  to  avoid  the  ICI,  a  COEDM  symbol  is  cyclically 
extended  by  adding  a  guard  time.  A  general  block  diagram  of  the  transmitter  and  the  re¬ 
ceiver  for  the  COEDM  scheme  is  shown  in  Eigure  3.  The  COEDM  is  used  in  802.1  la 
standard.  In  the  following  discussion  as  we  describe  the  COEDM  system,  frequent  refer¬ 
ences  will  be  made  to  this  standard. 
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Figure  3.  Block  Diagram  of  COFDM  Transmitter  and  Receiver  (From  Ref.  7.) 


1.  Channel  Coding 

In  the  IEEE  802.1  la  standard,  data  is  encoded  with  a  convolutional  encoder  with 
a  coding  rate  R  =  1/2,  2/3  or  3/4  .  The  convolutional  encoder  uses  the  industry-standard 
generator  polynomials,  gg  =  133g  and  gj  =  171g  of  rate  7?  =  1/2  [6].  The  code  rate  for  the 

convolutional  code  can  be  changed  by  using  a  puncturing  process.  In  this  thesis,  the  data 
are  also  encoded  with  RM  codes.  The  details  of  RM  codes  are  explained  in  Chapter  III. 


2.  Block  Interleaver 

If  decoding  errors  occur  in  a  codeword  and  these  are  passed  to  the  next  block, 
they  may  affect  the  performance  of  the  entire  system.  The  performance  of  the  system  can 
be  improved  if  these  errors  are  distributed  over  the  other  code  words.  This  can  be 
achieved  by  interleaver/deinterleaver. 

A  block  interleaver  consists  of  a  two-dimensional  array,  into  which  the  data  are 
read  along  its  rows.  When  the  array  is  full,  the  data  are  read  out  by  the  columns,  thus  the 
order  of  the  data  is  permuted.  The  original  order  can  be  received  by  the  corresponding 
deinterleaver  in  which  the  data  are  read  in  by  columns  and  read  out  by  rows. 
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3.  Symbol  Mapping 

The  interleaved  and  rearranged  data  are  mapped  onto  eonstellation  points  in  ac- 
cordanee  with  the  modulation  type.  Figure  4  shows  QPSK  constellation  points. 
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Figure  4.  QPSK  Constellation  points  (From  Ref.  7.) 

4,  Discrete  Fourier  Transform  (DFT) 

There  are  52  subcarriers  per  channel  in  a  IEEE  802.1  la  standard  WLAN  system, 
where  48  of  these  subcarriers  carry  data  and  the  remaining  four  subcarriers  are  used  as 
pilot  tones.  After  serial-to-parallel  conversion,  each  COEDM  symbol  is  modulated  over 
52  subcarriers  by  applying  an  inverse  fast  Eourier  transform  (lEET). 

5.  Guard  Interval  and  Cyclic  Extension 

A  guard  time  is  added  to  each  COEDM  symbol  to  eliminate  the  ISI  and  ICI,  and  it 
is  removed  before  the  LET  operation  at  the  receiver.  Since  the  other  parameters  are  cho¬ 
sen  according  to  the  guard  interval  time,  it  is  an  important  parameter  for  the  COEDM  sys¬ 
tem.  As  long  as  the  guard  time  is  larger  than  the  expected  delay  spread,  multipath  com¬ 
ponents  from  one  symbol  do  not  cause  interference  with  the  other  symbol. 
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A  zero-padding  extension  could  be  used  as  guard  interval;  however,  this  can  not 
eliminate  ICI  effects.  Therefore,  the  COFDM  symbol  is  cyclically  extended  into  the 
guard  interval.  The  cyclic  prefix  for  the  transmission  signal  consists  of  a  specific  part  of 
last  samples  that  are  added  to  the  beginning  of  the  signal.  Consequently,  we  have  an  ex¬ 
tended  signal.  To  ensure  orthogonality  between  different  subcarriers,  it  is  required  that  all 
subcarriers  differ  by  an  integer  number  of  cycles  within  the  FFT  integration  time.  Al¬ 
though  the  addition  of  guard  interval  eliminates  the  ICI  and  ISI  effects,  it  causes  ineffi¬ 
cient  use  of  the  frequency  spectrum 

6.  Modulator 

The  COFDM  signal  then  is  upconverted  to  the  5-GHz  band  and  transmitted  over 
the  channel. 


7.  Channel 

The  presence  of  noise  in  the  channel  affects  the  ability  to  make  correct  decisions 
about  the  received  symbols  at  the  receiver  part  of  the  communication  system,  thereby 
limiting  the  data  transmission  rate. 

For  additive  white  Gaussian  noise  (AWGN),  a  received  signal  is  typically  mod¬ 
elled  as  an  attenuated  desired  signal  in  AWGN  with  the  noise  having  a  power  spectral 
density  (PSD)  given  by 

G„(/)  =  ^W/Hz  (2-1) 

A  free  space  propagation  model  may  be  adequate  to  analyze  the  satellite  commu¬ 
nication  systems  since  there  is  always  a  line-of-sight  (LOS)  component  and  almost  no 
multipath  fading  effects.  However,  as  a  signal  travels  from  the  transmitter  to  the  receiver 
in  the  channels  of  interest,  the  path  it  takes  can  vary  from  a  simple  LOS  to  one  that  is  se¬ 
verely  obstructed  by  buildings,  hills,  and  trees,  which  lead  to  multipath  propagation.  The 
free  space  propagation  model  may  not  be  adequate  to  analyze  and  describe  the  character 
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istics  of  a  channel.  Because  of  the  uneven  terrain  and  other  obstructions,  a  multipath 
model  is  desired.  Consequently,  the  received  signal  undergoes  fading  as  it  may  consist  of 
multiple  waves,  which  have  randomly  distributed  amplitudes  and  phases. 

The  mean  excess  delay,  rms  delay  spread,  and  excess  delay  spread  are  parameters 
used  to  characterize  multipath  channels.  The  mean  excess  delay  (r  )  and  rms  delay 
spread  ( cr^ )  of  multipath  channels  are  used  to  quantity  the  time  dispersive  properties  of  a 
multipath  channel.  The  mean  excess  delay  is  the  first  moment  of  the  power  delay  profile 
as  given  by 


T  =  -  (2-2) 

k  k 

and  the  rms  delay  spread  is  the  square  root  of  the  second  central  moment  of  the  power 
delay  profile  defined  as 

(2-3) 

where 

=  - ,  (2-4) 

k  k 

where  is  absolute  power  level.  These  parameters  are  measured  relative  to  the  first  de¬ 
tectable  signal  arriving  at  Tq  =  0 .  Typical  values  of  rms  delay  spread  are  on  the  order  of 

microseconds  in  outdoor  channels  and  on  the  order  of  nanoseconds  in  indoor  channels. 

[8] 

Due  to  the  relative  motion  between  the  transmitter  and  the  receiver,  each  multi- 
path  signal  is  subjected  to  Doppler  shift.  The  Doppler  shift  is  computed  by 


10 


(2-5) 


V 

L=-fc 

c 

where  v  is  the  velocity  of  the  receiver,  c  is  the  speed  of  light  (c  =  3x10*  m/s),  and  is 
the  carrier  frequency.  [8] 

Different  transmitted  signals  undergo  different  types  of  fading  depending  on  sig¬ 
nal  and  channel  parameters.  Four  possible  effects  that  the  time  dispersion  and  frequency 
dispersion  in  a  channel  cause  are  frequency  selective  fading,  flat  fading,  slow  fading  and 
fast  fading. 


8,  Receiver 

On  the  receiver  side  of  the  COFDM  system,  the  reverse  operations  are  performed. 
At  the  front  end,  a  low-noise  amplifier  (LNA)  that  reduces  the  effective  noise  tempera¬ 
ture  of  the  receiver  and  an  automatic  gain  control  (AGC)  that  estimates  the  power  of  the 
pilot  tone  and  controls  the  power  at  the  demodulator  output  are  used.  The  guard  interval 
is  removed  once  the  symbols  are  detected.  The  symbol  constellations  are  recovered  by 
passing  the  signal  through  FFT.  The  resulting  data  are  deinterleaved  and  channel  de¬ 
coded. 


C.  PEAK  TO  AVERAGE  POWER  RATIO  (PAPR) 

COFDM  signals  are  generated  as  the  sums  of  a  number  of  subchannel  signals, 
which  are  continuous-time  sinusoidal  signals.  The  summation  may  cause  constructive 
interference,  thereby  resulting  in  a  high  peak-to-average  power  ratio  (PAPR). 

The  instantaneous  power  of  the  signal  is 

/>,(<)  =  |S,(0f  (2-6) 

where  S^(t)  is  the  transmitted  COFDM  signal. 

The  peak-to-average  power  ratio  (PAPR)  is  given  by 
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max|(5'^(0)|^ 

PAPR  = - ,  (2-7) 

^{|(*^.(0)|  } 

where  T  is  the  symbol  duration  and^{}  denotes  expectation. 

If  the  peak  transmission  power  is  physically  limited  or  limited  by  regulatory  rules, 
the  system  must  be  designed  with  reduced  average  transmitter  power.  This  may  reduce 
the  effective  transmission  range  of  the  system.  Additionally,  costs  may  rise  as  more 
equipment  may  be  needed  to  cover  the  same  transmission  range.  [9] 

In  brief,  high  PAPR  is  a  disadvantage  of  COFDM  communication  systems.  In  this 
thesis,  we  investigate  ways  to  reduce  PAPR. 


In  this  chapter,  the  COFDM  system  architecture  is  presented  and  the  advantages 
and  the  drawbacks  of  the  system  discussed.  The  PAPR  problem  in  COFDM  systems  is 
examined.  The  next  chapter  presents  the  coding  scheme  and  the  Reed-Muller  error  cor¬ 
rection  codes  as  a  mitigation  technique  to  reduce  the  required  PAPR. 
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III.  CODING  SCHEME 


This  chapter  begins  with  an  overview  of  error  correetion  eoding  in  digital  eom- 
munieation  systems  and  the  reason  why  ehannel  eoding  is  used.  Also,  we  diseuss  how 
error  correetion  eodes  improve  the  performanee  of  digital  eommunieation  systems  and 
present  a  summary  of  eommonly  used  error  eorreetion  eodes.  The  Reed-Muller  error  eor- 
reetion  eodes  are  diseussed  in  detail.  Also,  peak  elipping  and  Hanning  windowing  are 
explained. 


A,  ERROR  CORRECTION  SYSTEMS 

Due  to  the  effeets  of  noise  and  multipath  fading  in  the  ehannel,  the  transmitted 
signal  arrives  at  the  reeeiver  with  some  errors.  The  errors  in  the  demodulated  data  are 
eharaeterized  in  terms  of  a  BER,  whieh  is  direetly  proportional  to  the  symbol  rate  and 
inversely  proportional  to  transmitter  power  and  bit-energy  to  noise  power  speetral  density 
ratio  (Ej^INq  ).  The  bit  error  rate  is  an  important  performanee  parameter  of  digital  eom- 
munieation  systems. 

Error  eontrol  strategies  basieally  have  two  main  eategories,  automatie  repeat  re¬ 
quest  (ARQ)  and  forward  error  eorreetion  (EEC).  In  this  ehapter,  forward  error  eorreetion 
eoding  is  diseussed. 

In  forward  error  eorreetion  eoding,  a  eertain  number  of  redundant  bits  are  added 
to  data  bits  in  a  partieular  pattern  aeeording  to  the  type  of  the  eode.  In  other  words,  for 
every  k  data  bits,  n  eoded  bits  are  transmitted,  where  n>k  An  the  reeeiver,  the  k data 
bits  ean  be  reeovered  by  performing  a  deeoding  operation  on  the  n  reeeived  eoded  bits. 

The  transmission  eonditions  in  wireless  eommunieation  ehannels  are  severe  due 
to  multipath  fading  and  the  variation  of  the  signal-to-noise  power  ratio.  Therefore,  in  or¬ 
der  to  design  a  eommunieation  system  with  an  aeeeptable  BER,  error  eorreetion  eoding 
must  be  used  to  proteet  the  data  from  transmission  errors. 
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As  long  as  the  signal-to-noise  ratio  (SNR)  is  high  and  the  ehannel  is  relatively 
flat,  error  eorreetion  eoding  may  be  unneeessary  in  OFDM  systems.  However,  uneoded 
OFDM  systems  do  not  perform  well  in  fading  ehannels.  [10] 

The  code  rate  r  =  kjn  is  a  ratio  of  the  number  of  data  bits  to  the  total  number  of 
coded  bits  transmitted  per  code  word. 

As  a  result  of  redundancy,  the  number  of  bit  errors  can  be  expected  to  increase. 
However,  the  reliability  of  demodulated  data  is  increased  because  redundancy  is  used  to 
correct  some  of  the  errors.  A  better  BER  can  be  achieved  at  the  output  of  the  decoder  by 
using  an  appropriate  coding  scheme.  The  BER  improvement  provided  by  the  channel 
coding  is  generally  expressed  in  terms  of  the  required  to  achieve  the  same  per¬ 

formance  without  coding.  This  difference  in  E^lN^  is  called  coding  gain  as  given  by 


V  ^0  J uneoded.. 


V  "^0  J  coded.. 


(3-1) 


Using  commercially  available  standard  error  correction  systems,  coding  gains  up 
to  6  to  9  dB  are  achievable.  [11,  12] 

Einear  block  codes,  convolutional  codes  and  turbo  codes  are  the  most  widely  used 
error  correction  codes.  Although  only  convolutional  codes  and  Reed-Muller  codes  are 
used  in  the  COEDM  simulation  undertaken  in  this  thesis,  linear  block  codes  and  turbo 
codes  are  also  mentioned. 


Einear  block  codes  generate  n  coded  bits  for  k  data  bits  where  n>k  .  Eirst,  k  data 
bits  are  transmitted  and  then  n-k  redundant  bits  which  are  generated  according  to  the 
generator  matrices  of  the  code  are  transmitted.  The  encoder  consists  of  a  k-stage  shift  reg¬ 
ister  and  n-k  modulo-2  adders.  The  shift  register  outputs  are  connected  to  modulo-2 
adders  according  to  the  generator  matrices  of  the  code. 

The  main  difference  between  block  codes  and  convolutional  codes  is  that  a  sys¬ 
tem  utilizing  block  codes  transmits  the  k  data  bits  unaltered  and  then  transmits  the  n-k 
redundant  bits.  A  system  utilizing  convolutional  codes  produces  n  coded  bits  from  k 
data  bits,  and  the  code  word  does  not  contain  unaltered  k  data  bits. 
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As  mentioned  in  Chapter  II,  in  the  IEEE  802.1  la  WEAN  system,  the  data  are 
eoded  using  a  eonvolutional  eneoder  with  eoding  rate  r  =  1/2,  2/3  ,  or  3/4  .  The  eonvolu- 
tional  eneoder  uses  the  half-rate  industrial-standard  generator  polynomials,  =  133g  and 
gj  =  IVlg .  Higher  rates  are  obtained  from  the  half-rate  eode  by  using  a  puneturing  proe- 

ess,  whieh  omits  some  of  the  eneoded  bits  in  the  transmitter  and  inserts  a  dummy  zero  in 
plaee  of  omitted  bits  in  the  eonvolutional  deeoder  in  the  reeeiver.  [7] 

A  eoneatenated  eode  eonsists  of  two  separate  eodes  in  series  in  whieh  the  first 
eode,  ealled  outer  eode,  direetly  takes  the  information  bits  and  encodes  them.  The  second 
code,  called  inner  code,  takes  the  bits  coded  by  the  outer  code  and  encodes  them.  In  order 
to  prevent  decoding  errors  from  passing  from  one  coder  to  the  other,  the  errors  are  dis¬ 
tributed  by  using  the  interleaving/deinterleaving  operation.  The  interleaver  is  placed  be¬ 
tween  the  outer  and  inner  encoders  of  a  concatenated  code  in  the  transmitter  and  the  dein¬ 
terleaver  is  placed  between  the  inner  and  outer  decoders  in  the  receiver.  To  make  the  cod¬ 
ing  system  more  efficient,  the  output  of  the  outer  decoder  is  reapplied  to  the  inner  de¬ 
coder.  This  is  the  basis  of  the  iterative  decoding.  It  can  be  summarized  as  reapplying  the 
decoded  word  not  just  to  the  inner  code,  but  also  to  the  outer  and  repeating  as  many  times 
as  necessary  [13].  Turbo  codes  are  based  on  parallel  convolutional  concatenated  codes 
(PCCC).  Additionally,  pseudorandom  interleavers  instead  of  rectangular  interleavers  are 
used. 

Since  the  goal  of  this  thesis  was  to  analyze  the  COEDM  system  with  Reed-Muller 
(RM)  error  correction  codes,  the  focus  was  on  the  encoding  and  decoding  algorithms  of 
RM  codes. 

B.  REED-MULLER  CODES 

This  section  describes  the  encoding/decoding  algorithm  of  the  Reed-Muller  (RM) 
coding  and  how  it  is  used  for  the  reduction  of  PAPR  in  COEDM  systems.  Since  the  en¬ 
coding  and  decoding  algorithms  are  complicated  and  different  from  the  other  schemes, 
some  binary  operations  used  with  RM  codes  are  first  defined  and  then  the  encoding  and 
decoding  algorithms  are  presented. 


15 


RM  codes  were  described  by  Muller  in  1954,  and  Reed  provided  a  better  alge- 
braie  representation  with  a  deeoding  algorithm  in  the  same  year.  RM  eodes  were  used  in 
spaee  applieations  between  1969  and  1977.  Although  they  seem  to  have  lost  their  attrae- 
tion  in  the  spaee  program  beeause  of  the  adoption  of  the  newer  eodes,  it  does  not  neees- 
sarily  mean  that  BCH  (Bose-Chaudhuri-Hoequenghem)  and  Reed-Solomon  eodes  are  a 
better  choice  relative  to  RM  eodes  for  all  applieations.  [11] 

1.  Definitions  of  Operations 

Three  possible  operations  performed  in  encoding  and  decoding  are  addition,  mul¬ 
tiplication  and  dot  product.  For  the  binary  vectors  x  =  . ,x^)  and 

y  =  (yi,y2’ . ’T«)  >  where  each  x.  or  is  either  1  or  0,  addition  is  defined  as  [1] 

x+y={x,+y„X2+y2, . (3-2) 

and  the  complement  x  of  veetor  x  is  the  vector  equal  to  1  +  x .  Multiplieation  is  defined 
as 

X*  y={x,^y„X2*y2, . ,x„*y„).  (3-3) 

Finally,  the  dot  product  of  x  and  y  is  defined  as 

x.y=(Xi*yi+X2*y2+ . +  (3-4) 

2.  Encoding  Algorithm 

The  RM  coding  scheme  is  denoted  in  terms  of  R{r,m) .  The  parameters,  r  and  m, 
are  defined  by 

m  =  log2  n  (3-5) 

where  n  is  the  number  of  the  columns  of  the  eneoding  matrix  and 

r  =  m-log2t/„i„  (3-6) 

where  is  the  minimum  Hamming  distance.  The  number  of  the  rows  of  the  eneoding 
matrix  is  given  by  [1 1] 
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(3-7) 


where  x-  is  the  monomial  vector  and  has  a  pattern  of  2“  '  ones  followed  by  I""  ‘  zeros, 
repeated  until  n  is  completed.  The  first  row  of  the  encoding  matrix  is  all  ones  or  all 

zeros  and  the  sequence  of  ones  or  zeros  in  the  next  rows  are  interchangeable.  In  this  the¬ 
sis,  the  first  row  is  all  ones  and  the  sequence  pattern  is  ones  and  zeros,  respectively.  For 
example,  if  m  =  3,  then  the  encoding  matrix  for  R(l,3)  is  given  by 

1  1  1  1  1  1  1  11  1 

1  1  1  1  0  0  0  0  x, 

(3-9) 

1  1  0  0  1  1  0  0  X2 

1  0  1  0  1  0  1  oJx3 

where  the  column  on  the  right  hand  side  is  the  corresponding  monomial  vector. 
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The  R(2,3)  code  can  be  generated  by  adding  the  rows  XjXj ,  XjX3  and  x^x^  to  the 


generator  matrix  in  (3-9)  as  follows: 


G(2,3) 


1 

1 

1 

1 

1 

1 

1 

r 

1 

1 

1 

1 

1 

0 

0 

0 

0 

Xi 

1 

1 

0 

0 

1 

1 

0 

0 

X2 

1 

0 

1 

0 

1 

0 

1 

0 

X3 

1 

1 

0 

0 

0 

0 

0 

0 

X1X2 

1 

0 

1 

0 

0 

0 

0 

0 

XjX3 

1 

0 

0 

0 

1 

0 

0 

0) 

X2X3 

(3-10) 


where  the  column  in  the  right  hand  side  is  the  corresponding  monomial  vector.  The  gen¬ 


erator  matrix  for  R(3,3)  is  obtained  by  adding  the  x^X2XJ  term  to  the  above  and  is  given 
by 


G(3,3) 


1  1 
1  1 
1  1 
1  0 
1  1 
1  0 
1  0 
1  0 


1  1 
1  1 
0  0 
1  0 
0  0 
1  0 
0  0 
0  0 


1  1 
0  0 
1  1 
1  0 
0  0 
0  0 
1  0 
0  0 


1  1 
0  0 
0  0 
1  0 
0  0 
0  0 
0  0 
0  0 


1 

Xi 

X2 

X3 

XjXj 

XjX3 

X2X3 

XjX2X3 


(3-11) 


It  must  be  noted  that  the  R{r,m)  encoding  matrix  can  be  achieved  by  adding 


rows  to  the  R{r-\,m)  encoding  matrix.  Since  the  encoding  matrix  has  k  rows,  the 


message  is  sent  to  the  encoding  matrix  in  blocks  of  length  k  .  If  the  message  is  defined  as 


m=  . 


(3-12) 


then  the  encoded  message  is 

i=l 

where  R.  is  the  corresponding  row  of  the  encoding  matrix. 


(3-13) 
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3,  Decoding  Algorithm 

The  decoding  algorithm  is  somewhat  complicated  and  consists  of  a  few  steps.  It  is 
based  on  majority  logic  or  threshold  decoding.  In  general,  majority  logic  techniques  are 
fast  and  easy  to  implement. 

An  R{r,  m)  code  can  correct  up  to  e  errors  if  the  distance  between  two  code 
words  is  greater  than  2e,  where  the  distance  is  described  as  the  number  of  different  values 
at  the  corresponding  places  between  two  code  words. 

Characteristic  vectors  are  the  vectors  whose  dot  product  is  one  with  the  corre¬ 
sponding  row  and  zero  with  all  the  other  rows  of  the  encoding  matrix.  The  decoding  algo¬ 
rithm  begins  with  finding  2"'^''  characteristic  vectors  for  each  row  of  the  encoding  matrix. 
After  finding  the  characteristic  vectors,  the  next  step  is  the  dot  product  of  these  character¬ 
istic  vectors  with  the  received  encoded  message  .  This  operation  gives  2“ ''  results 

and,  by  performing  the  majority  logic,  the  coefficient  of  the  row  is  calculated.  In  other 
words,  if  the  number  of  the  zeros  is  larger  than  the  number  of  ones  in  the  dot  product  op¬ 
erations,  the  coefficient  is  zero  and  one  otherwise.  These  steps  thus  far  are  performed  for 
every  row  of  the  encoding  matrix  except  the  first  one.  The  coefficients  are  then  multi¬ 
plied  with  the  corresponding  row  and  the  results  are  added.  At  this  point,  is  obtained, 

which  is  a  row  vector  having  2"'  columns.  The  addition  of  and  provides  the  loca¬ 
tion  of  the  error(s)  and  knowledge  about  the  first  row  of  the  encoding  matrix,  that  is, 
whether  it  is  all  ones  or  all  zeros.  If  the  number  of  the  ones  in  this  row  is  larger  than  the 
number  of  zeros,  the  first  row  of  the  matrix  is  all  ones,  otherwise  it  is  all  zeros.  The  1  ’s 
and  O’s  so  obtained  form  the  recovered  message.  [1] 

It  must  be  noted  that  the  error-correcting  capability  of  the  RM  codes  can  be  in¬ 
creased  by  increasing  the  minimum  Hamming  distance. 

4.  Reduction  of  PAPR  of  COFDM  Using  RM  Coding 

While  discussing  the  drawbacks  of  the  COFDM  system,  it  was  mentioned  that 
high  PAPR  is  the  main  disadvantage  of  the  system. 
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Some  research  has  been  performed  in  order  to  reduce  the  PAPR  in  the  COFDM 
systems.  One  of  the  solutions  to  reduce  PAPR  requires  the  use  of  block  codes.  This  ap¬ 
proach  is  based  on  the  selection  of  appropriate  code  words  that  reduce  the  PAPR,  and 
then  the  use  of  these  code  words.  However,  performing  a  long  search  to  find  these  appro¬ 
priate  codes  and  the  necessity  for  large  lookup  tables  make  this  approach  inefficient  [14]. 
Another  approach  is  based  on  introducing  a  specific  phase,  which  is  independent  and 
known  to  both  the  transmitter  and  the  receiver  for  each  coordinate  of  all  the  code  words. 
This  approach  obtained  4.5  dB  in  PAPR  reduction  by  using  the  computed  phase  shifts 
[10].  In  this  work,  the  use  of  RM  codes  in  order  to  reduce  PAPR  is  explored. 

Davis  and  Jedwab  in  [15]  showed  that  large  sets  of  binary  length  of  2“  Golay 
complementary  pairs  can  be  obtained  from  the  certain  second-order  cosets  of  the  first  or¬ 
der  RM  code.  These  cosets  can  be  grouped  according  to  their  PAPR  values  for  the  codes 
of  lengths  I"" ,  where  m  =  2,3,4,5  .  The  first  group  of  cosets  has  a  PAPR  of  at  most  2. 

The  second  one  has  a  PAPR  of  at  most  4  and  so  forth.  An  effective  combination  of  the 
code  rate,  PAPR  and  minimum  distance  can  be  obtained  by  selecting  second-order  cosets 
from  the  list.  This  method  can  be  efficiently  used  for  the  COFDM  system  for  a  small 
number  of  carriers.  For  a  larger  number  of  carriers,  it  becomes  inefficient  since  the 
evaluation  of  PAPR  values  of  a  large  set  of  cosets  requires  a  large  amount  of  computa¬ 
tion.  [15] 

In  summary,  using  RM  codes  with  the  COFDM  system  is  expected  to  not  only 
provide  efficient  encoding  and  decoding  algorithms  and  high  error  correction  capability 
with  high  code  rates,  but  also  to  reduce  PAPR,  especially  for  a  small  number  of  carriers. 


C.  PEAK  CLIPPING 

Peak  clipping  is  a  simple  approach  to  reduce  the  PAPR.  Since  the  peaks  in  the 
COFDM  signal  occur  with  a  low  probability,  peak  clipping  can  be  considered  an  effec¬ 
tive  technique  for  the  reduction  of  PAPR.  However,  it  must  be  noted  that  it  is  a  nonlinear 
process  and  may  worsen  the  BER  performance  due  to  the  inband  signal  distortion.  Addi¬ 
tionally,  it  increases  the  out-of-band  radiation  and,  therefore,  reduces  the  frequency  spec¬ 
trum  efficiency.  The  clipping  ratio,  CR,  is  defined  by 
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(3-14) 


cr 


where  A  is  the  elip  level  and  cr  is  the  rms  power  of  the  COFDM  signal.  The  COFDM 
signal  is  clipped  whenever  it  exceeds  the  clip  level. 


D.  HANNING  WINDOWING 

Hanning  windowing  is  another  approach  to  reduce  the  PAPR  of  COFDM  signals. 
The  COFDM  signal  is  multiplied  by  the  window  function  when  the  signal  exceeds  the 
clipping  level  or  falls  below  the  bottom  level.  While  the  clipping  operation  directly  chops 
off  the  peaks,  windowing  results  in  a  smooth  signal.  The  peak  window  is  given  by 

ri-k^(0.5-0.5cos(2;rn/M))  0<n<M 
wSn]  =  \  .  (3-15) 

0  otherwise 


where  is  the  attenuation  factor  of  the  window  and  M  is  the  width  of  the  window.  The 
bottom  window  is  given  by 


Jl  +  (0.5  -  0.5  cos(2;rn  /  M)) 

[o 


0<n<M 

otherwise 


(3-16) 


where  is  the  amplification  factor  of  the  window.  Figure  9  shows  a  9-point  peak  and  a 
9-point  bottom  Hanning  windows.  [2] 


9*point  Peak  Hanning  Window 


9-point  Bottom  Hanning  Window 


Figure  5.  9-Point  Hanning  Windows 
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When  the  windows  are  multiplied  with  the  COFDM  signal,  the  resulting  speetrum 
is  the  speetrum  of  the  windowed  signal.  Therefore,  the  window  width  is  an  important  pa¬ 
rameter  that  affeets  the  BER  performaneeof  the  system. 


This  ehapter  discussed  the  error  correcting  coding  schemes.  The  RM  coding 
scheme  was  presented  in  detail.  The  RM  coding  is  known  to  reduce  the  PAPR,  which  is  a 
problem  of  the  COFDM  systems.  The  encoding  and  decoding  algorithms  of  RM  codes 
were  presented.  Peak  clipping  and  Hanning  window  techniques  were  explained.  The  next 
chapter  presents  the  COFDM  simulation  architecture  and  reports  the  results  of  simula¬ 
tion. 
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IV.  SIMULATION  RESULTS 


This  chapter  describes  the  simulation  model  used  in  this  thesis  and  presents  the 
simulation  results.  The  MATLAB  simulation  code  in  [1]  was  modified  and  expanded  to 
implement  the  simulation  model.  Simulations  were  performed  to  study  the  system  per¬ 
formance  with  RM  codes  under  a  variety  of  channel  conditions.  The  system  was  also 
tested  for  PAPR  reduction.  In  addition  to  the  RM  codes,  peak  clipping  and  Hanning  win¬ 
dowing  were  implemented  in  order  to  reduce  the  PAPR. 


A,  SIMULATION  MODEL 

Since  the  main  goal  of  this  thesis  was  to  simulate  the  COFDM  system  by  utilizing 
RM  codes,  all  signal  processing  including  forward  error  correction  coding  (FEC)  and  dif¬ 
ferent  types  of  channel  characteristics  are  performed  at  the  base  band.  The  block  diagram 
of  the  entire  system  is  shown  in  Figure  6.  A  detailed  description  of  the  different 
MATLAB  functions  corresponding  each  of  the  blocks  in  Figure  6  is  provided  in  Appen¬ 
dix  A. 


Figure  6.  COFDM  System  Block  Diagram. 


B,  SIMULATION  PARAMETERS 

During  the  simulations,  in  order  to  compare  the  results,  the  same  random  message 
was  used.  The  seed  value  in  msg.m  was  changed  only  for  comparing  the  PAPR  values  of 
the  system  with  different  random  messages.  Although  the  code  makes  it  possible  to  select 
a  different  number  of  symbols  and  interleaver  pairs,  all  simulation  runs  were  performed 
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with  12,000  symbols  and  a  (120,  100)  interleaver  pair.  While  testing  the  eonvolutional 
eodes,  the  IEEE  802.1  la  standard  half-rate  eonvolutional  eodes  with  A:  =  50  and  n  =  100 
were  performed. 

After  eonstruetion  of  the  subbloeks  and  modifieation  of  the  eode  in  [1],  five  dif¬ 
ferent  ehannels  were  formed.  Channel  0  is  a  noise-free  ehannel  with  no  AWGN  and  mul¬ 
tipath  effeets.  Channels  1,  2  and  3  ineorporate  AWGN,  multipath,  and  AWGN  and  multi- 
path,  respeetively.  Channel  4  ineorporates  the  outdoor  ehannel  eharaeteristies  of  AWGN, 
severe  multipath,  and  mobility. 

In  Channel  1 ,  the  standard  deviation  <j  of  white  Gaussian  noise  is  varied  from  0 
to  0.06  for  different  eoding  options. 

The  multipath  fading  parameters  used  in  Channels  2  and  3  are  tabulated  in  Table 
1 .  The  multipath  loss  in  dB  and  the  delay  in  ms  are  listed  for  indoor  ehannel  environment. 
There  are  18  taps  and  delay  eoeffieients  in  the  ehannel.  Three  different  Doppler  frequen- 
eies  of  5  Hz,  10  Hz,  and  15  Hz  were  eonsidered;  the  eorresponding  veloeities  of  these 
frequeneies  were  0.29,  0.58  and  0.87  m/s,  respeetively.  They  represent  walking  speeds  in 
an  indoor  environment. 


Loss  (dB) 

[0,  2.17,  4.34,  6.51,  8.69,  10.86,  13.03,  15.20,  17.37, 
19.54,  21.71,  23.89,  26.06,  28.23,  30.4,  32.57,  34.74, 
36.92] 

Delay  (msec) 

[0,  0.05,  0.10,  0.15,  0.20,  0.25,  0.30,  0.35,  0.40,  0.45, 

0.50,  0.55,  0.60,  0.65,  0.70,  0.75,  0.80,  0.85] 

Doppler  (Hz) 

5 

10 

15 

Table  1 .  Multipath  Eading  Parameters  for  Channels  2  and  3 


The  outdoor  ehannel  parameters  of  Channel  4  are  shown  in  Table  2.  Within  this 
ehannel,  two  different  eases  eaeh  with  different  loss,  delay  and  Rieian  faetors  are  imple¬ 
mented. 


Channel  4A 

Channel  4B 

Loss  (dB) 

0.0,  1.0,  9.0,  10.0,  15.0, 
20.0 

2.5,  0.0,  12.8,  10,  25.2,  16 

Delay  (msec) 

0.0,0.25,0.5,  1.0,  1.9, 

2.2 

0.0,0.25,9.0,  13.0,  17.0,  20.0 

Rician  Factors  (%) 

0,  0,  0,  0,  0,  0 

0.5,  0.5,  0,  0,  0,  0 

Table  2.  Outdoor  Channel  Parameters  of  Channel  4 
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In  order  to  analyze  the  effects  of  using  RM  codes  in  the  COFDM  system  for  both 
indoor  and  outdoor  channel  characteristics,  a  step-by-step  simulation  methodology  was 
followed.  The  performance  of  the  COFDM  system  with  RM  codes  is  tested  under  differ¬ 
ent  conditions  in  the  order  presented  in  Table  3. 


Test  No, 

Channel 

Channel  Description 

1 

Channel  0 

Noise-free  Channel 

2 

Channel  1 

AWGN  effect 

3 

Channel  2 

Multipath  effect 

4 

Channel  3 

AWGN  +  Multipath  effect 

5 

Channel  4 

Severe  mobile  outdoor  channel 

Tables.  General  Test  Plan. 


C.  SIMULATION  RESULTS 

Channel  0  was  used  to  test  whether  the  system  was  configured  properly  and  work¬ 
ing  correctly.  In  this  case,  the  received  message  is  the  same  as  the  source  message.  Nu¬ 
merous  simulations  performed  for  different  types  of  RM  and  convolutional  codes  demon¬ 
strated  that  the  code  ran  correctly.  Figure  7  shows  the  transmitted  and  received  QPSK 
constellations. 
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0.5 


-0.5 


-I.5L- 

-1.5 


■  1  I 


-1  -0.5 


0.5  1  1.5 


t 

- H 

— 

- - 

— 

— 

_ 

.5 1 - 1 - 1 - 1 - 1 - 1 - 

-1.5  -1  -0.5  0  0.5  1  1.5 


(a)  Transmitted 


(b)  Received 


Figure  7.  Transmitted  and  Received  QPSK  Signal  Constellations  for  Channel  0. 
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The  performance  of  the  COFDM  system  is  next  tested  using  Channel  1 ,  which  in¬ 
cludes  AWGN  without  multipath  fading  effects.  The  received  QPSK  constellations  for 
different  levels  of  noise  are  shown  in  Figure  8. 
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1 

► 

. H 

. 

. 

i 

*• 

. 

5  I - i - 1 - 1 - i - 1 - 

-1.5  -1  -0.5  0  0.5  1  1.5 

0-  =  0,009 


a  =  0.02 


Figure  8.  The  Effects  of  AWGN  over  QPSK  Signal  Constellation. 


The  BER  plots  were  obtained  for  different  types  of  RM  codes  and  convolutional 
codes.  The  number  of  errors  that  an  R{r,  m)  code  can  correct  depends  on  the  r  and  m 
values.  The  number  of  the  errors  that  different  types  of  RM  codes  can  correct  is  tabulated 
in  Table  4.  The  bit-error  plots  for  three  different  RM  codes  and  a  rate  1/2  convolutional 
code  along  with  the  QPSK  theoretical  plots  are  shown  in  Eigure  9. 
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Type  of 
code 

Number  of  the  rows 
(number  of  the  bits  in  a 
block) 

Number  of  the  er- 
ror(s)  that  can  be 
corrected 

Rate 

R{r,m) 

( m\  ( m\  ( m\ 

—  1+  +  +....+ 

U  y  v2  y  y 

2'”^''^i-l 

2m-,-l_i 

k 

R(l,3) 

4 

1 

-  =  0.25 

4 

R(2,4) 

11 

1 

—  =  0.0909 

11 

R(2,5) 

16 

3 

—  =  0.1875 

16 

able  4.  The  Number  of  the  Error(s)  that  Different  Types  of  RM  Codes  Can  Correct. 


1 Q-'* - ^ ^ - 1 - 1 - ^ - 

2  4  6  8  10  12  14 

Eb/No(dB) 

Figure  9.  BER  Performance  of  the  COFDM  system  under  Channel  1  (AWGN  only) 

Conditions. 


From  Figure  9  and  Table  4,  the  performance  of  an  RM  code  depends  on  the  num¬ 
ber  of  error(s)  that  it  can  correct  in  a  block.  Since  this  rate  is  worst  for  R  (2,4),  the  corre¬ 
sponding  BER  performance  curve  was  also  the  worst.  Compared  to  the  theoretical,  there 
is  2.5  dB  difference  for  convolutional  code,  2.7  dB  for  R(l,3),  3.6  dB  for  R(2,5)  and  4.8 
dB  for  R(2,4)  at  =  10  ^ 
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1 .  Multipath  Effects 

The  performance  of  the  COFDM  system  was  studied  by  adding  multipath  effects 
to  the  channel.  First,  multipath  effects  without  AWGN  were  simulated  using  Channel  2. 
Figure  10a  shows  the  received  QPSK  constellation.  The  constellation  points  are  scattered 
from  their  original  position  due  to  the  effects  of  multipath  fading.  Figure  10b  shows  that, 
with  the  addition  of  differential  decoding,  the  constellation  points  realigned  somewhat 
within  their  respective  spaces. 


20 
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5 

0 

-5 
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-15 

-20 

-5  0  5  -20  -15  -10  -5  0  5  10  15  20 

(a)  Before  differential  decoding  (b)  After  differential  decoding 

Figure  10.  The  Effects  of  Multipath  on  QPSK  Signal  Constellation 


Channel  3  takes  AWGN,  multipath  effects,  and  Doppler  shift  into  account,  hence 
more  realistic  than  Channel  2.  Channel  3  may  be  considered  a  good  representation  for 
indoor  environments.  Figure  1 1  shows  the  magnitude  of  the  QPSK  signal  at  the  input  of 
the  receiver.  The  magnitude  is  plotted  for  different  symbols  at  different  subcarriers  for  an 
AWGN  standard  deviation  of  cr  =  0.02  .  Note  that  ideally  the  magnitude  is  a  constant  for 
all  symbols  at  all  subcarriers.  Figures  12  and  13  show  the  QPSK  constellation  prior  to 
and  after  the  differential  decoding.  QPSK  constellation  points  are  shifted  from  their 
original  phase  sectors  because  of  the  effects  of  multipath  fading.  The  spreading  of  the 
constellation  points  increases  with  the  noise  variance.  As  Figure  13  illustrates,  the  con¬ 
stellation  points  can  be  realigned  to  their  respective  phase  sectors  to  some  extent  by  using 
differential  decoding. 
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Symbol  Row  Number  0  OFDM  Freq  H 

a  =  0.02 

Figure  1 1 .  The  Representation  of  COFDM  Signal  Under  Multipath  and  AWGN  Ef¬ 
fects. 
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Figure  12.  The  Effects  of  AWGN  and  Multipath  on  QPSK  Signal  Constellation  Be¬ 
fore  Differential  Decoding 
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<7  =  0.001 


<7  =  0.02 


Figure  13.  The  Effeets  of  AWGN  and  Multipath  on  QPSK  Signal  Constellation  After 

Differential  Decoding 


Figures  14-16  show  the  bit-error  rate  plots  for  three  different  Doppler  shift  values. 
Each  figure  contains  three  different  types  of  RM  codes  and  a  rate  1/2  convolutional  code. 
The  BER  performance  curves  show  that  the  effect  of  the  5 -Hz  Doppler  frequency  creates 
slight  difference  in  the  BER  performance  of  the  system.  As  seen,  the  preformence  of  the 
convolutional  coder  is  better  than  that  of  the  RM  codes. 
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Figure  14.  BER  Performance  of  the  COFDM  system  under  Channel  3  (AWGN  + 

Multipath  Effect)  with  5-Hz  Doppler  Conditions 


Figure  15.  BER  Performance  of  the  COFDM  system  under  Channel  3  (AWGN  + 

Multipath  Effect)  with  10-Flz  Doppler  Conditions 
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Figure  16.  BER  Performance  of  the  COFDM  system  under  Channel  3  (AWGN  + 

Multipath  Effect)  with  15-Hz  Doppler  Conditions 


2,  Outdoor  Environment 

We  now  study  the  COEDM  system  performance  under  mobile  outdoor  channel 
conditions.  Two  channels  were  considered  for  this  purpose:  Channels  4A  and  4B. 

Figure  17  shows  the  COFDM  signal  magnitude  for  Channels  4A  and  4B  for 
<j  =  0.02  .  Although  the  transmitted  COFDM  signal  magnitude  was  at  unity,  the  multi- 
path  fading,  delay,  and  Rician  factors  in  the  channels  resulted  the  noticeable  signal  varia¬ 
tions  in  received  signal  power  levels.  The  corresponding  QPSK  constellations  are  shown 
in  Figure  18.  Figure  18  shows  that  constellation  points  were  scattered  from  their  original 
points  due  to  the  multipath  fading,  delay,  and  Rician  factors  of  the  channels.  Figure  19 
shows  the  QPSK  constellations  after  the  differential  decoding  for  Channels  4A  and  4B. 
Compared  to  the  constellation  points  in  Figure  12,  which  includes  multipath  effects  in 
AWGN,  the  signal  constellation  points  in  Figure  18  underwent  more  severe  multipath 
fading. 
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Symbol  Row  Number 


0  0 


OFDM  Freq  # 


Symbol  Row  Number 


0  0 


OFDM  Freq  # 

(a)  Channel  Model  4A  (b)  Channel  Model  4B 

Figure  17.  The  Effeets  of  Channel  4  on  Signal  Magnitude. 
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(b)  Channel  Model  4B 

Signal  Constellation  (Before  Differential  De¬ 
coding). 


ire  18. 


(a)  Channel  Model  4A 

The  Effects  of  Channel 
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(a)  Channel  Model  4A 


(b)  Channel  Model  4B 


Figure  19.  The  Effects  of  Channel  4  over  Signal  Constellation  (After  Differential  De¬ 
coding). 

Figures  20  and  21  show  the  BER  plots  for  Channel  4.  In  order  to  obtain  a  BER 
less  than  10'^,  the  value  should  be  higher  than  about  10  dB  for  RM  codes  and  5 

dB  for  convolutional  codes.  This  result  shows  that  the  COFDM  system  requires  high  sig- 
nal-to-noise  ratio  in  the  outdoor  environments. The  necessity  of  high  transmission  power 
for  Channels  4A  and  4B  creates  the  PAPR  problem.  The  results  of  the  PAPR  values  for 
the  system  are  presented  in  the  next  section. 
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BIT  Error  Rate 


Figure  20.  BER  Performance  Curves  for  Channel  4A  Simulation. 


Eb/No(dB) 


Figure  21 .  BER  Performance  Curves  for  Chaimel  4B  Simulation. 

35 


3.  PAPR  Reduction 

The  PAPR  values  of  the  COFDM  system  simulated  using  different  seeds  to  gen¬ 
erate  the  random  messages  are  tabulated  in  Table  5.  The  RM  eodes  used  in  this  work  are 
R(l,3),  R(2,4),  and  R(2,5).  The  RM  eodes  in  general  are  not  eonsistent  in  redueing  the 


PAPR  in  COFDM.  This  result  does  not  agree  with  the  elaim  in  [15]. 


FEC 

Code 

Type 

PAPR  (dB) 

Seed  =  13 

Seed  =  23 

Seed  =  33 

Average 

Max 

Ave 

Max 

Ave 

Max 

Ave 

Max 

Ave 

R(l,3) 

9.4734 

6.7845 

10.3031 

6.6575 

10.7564 

9.9686 

10.1776 

7.8035 

R(2,4) 

9.2566 

6.7392 

9.5706 

6.9064 

9.5816 

6.8251 

9.4696 

6.8235 

R(2,5) 

11.4281 

6.7306 

9.6601 

6.7231 

9.9839 

6.6464 

10.3573 

6.7000 

CONV. 

10.0090 

6.9200 

9.4433 

6.5059 

10.0320 

6.7182 

9.8281 

6.5347 

Table  5 .  PAPR  Values  without  Windowing  and  Peak  Clipping. 


Peak  elipping  and  Hanning  windowing  were  used  to  reduee  the  PAPR.  It  ean  be 
seen  that  peak  elipping  redueed  PAPR  by  3  dB.  However,  the  redueed  PAPR  value  intro- 
dueed  a  high  BER  as  ean  be  seen  in  Figure  22.  This  presents  a  tradeoff  between  PAPR 
reduction  and  BER  performance. 

Table  6  shows  the  results  when  the  COEDM  signal  was  clipped  at  18  through  12. 
Erom  Eigure  22,  the  BER  performance  of  the  system  at  CE  =  18  is  very  close  to  no¬ 
clipping  performance.  However,  as  the  peak  clipping  level  is  decreased,  the  required 
Ej  jN^  to  achieve  the  same  BER  performance  increased  because  of  the  increasing  prob¬ 
ability  of  existence  of  the  COEDM  signal  magnitudes  higher  than  the  clipping  level. 


PAPR(dB) 

Seed=33 

Max 

average 

No  Clipping 

10.7564 

9.9686 

CL=18 

10.4895 

6.9636 

CL=16 

9.5959 

6.9440 

CL=14 

8.5536 

6.8727 

CL=12 

7.3191 

6.4491 

Table  6.  Peak  Clipping  Results. 
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Figure  22.  BER  Performance  of  the  System  with  Peak  Clipping. 


Hanning  windowing  was  also  implemented  to  reduce  the  PAPR.  The  values  of 
=  0.1  and  =  0.2  were  assigned  as  Hanning  windows  coeffiecients  during  the  simu¬ 
lations  to  achieve  the  best  performance. 

When  the  windows  were  applied  to  the  COFDM  signal,  the  resulting  spectrum 
was  the  spectrum  of  the  windowed  signal.  The  windowing  process  was  tested  with  3-,  5- 
and  9-point  Hanning  windows.  The  PAPR  values  are  shown  in  Table  7  and  BER  per¬ 
formance  plots  are  shown  in  Figures  23  through  25  for  3-,  5-,  and  9-point  windows,  re¬ 
spectively.  Table  7  shows  that  the  improvement  in  PAPR  reduction  is  limited  with  Han¬ 
ning  windowing.  Figures  23  through  25  show  that  as  the  width  of  the  window  is  in¬ 
creased,  the  bit  errors  increased  as  well. 
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PAPR(dB) 

Seed  =  33 

Seed  =  23 

3-point  Hanning 
Windowing 

5-point  Hanning 
Windowing 

9-point  Hanning 
Windowing 

9-point  Hanning 
windowing 

Max 

Ave 

Max 

Ave 

Max 

Ave 

Max 

Ave 

No  win¬ 
dowing 

10.7564 

6.9686 

10.7564 

6.9686 

10.7564 

6.9686 

10.3031 

6.6575 

CR=1,4 

9.9869 

6.4219 

10.0097 

6.5186 

10.0486 

6.4939 

9.6236 

6.6020 

CR=1,4 

&B=1 

9.9144 

6.4307 

9.8565 

6.4740 

9.7377 

6.5390 

9.1957 

6.2795 

CR=1,4 
&  B=1.5 

9.9004 

6.4102 

9.8284 

6.4701 

9.6647 

6.5507 

11.4741 

6.3601 

CR=1,4 

&B=2 

9.8716 

6.4346 

9.7694 

6.5206 

9.7151 

6.6559 

11.3221 

6.4451 

Table  7.  The  Results  of  3,  5,  9-Point  Hanning  Windowing. 


Eb/No(dB) 

Figure  23.  BER  Performance  of  the  System  with  3-Point  Hanning  Windowing  (CR; 

Clipping  Ratio,  B;  Bottom  Level) 
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Figure  24.  BER  Performance  of  the  System  with  5-Point  Hanning  Windowing  (CR; 

Clipping  Ratio,  B:  Bottom  Level) 


Figure  25.  BER  Performance  of  the  System  with  9-Point  Hanning  Windowing  (CR; 

Clipping  Ratio,  B;  Bottom  level) 
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This  chapter  studied  the  simulation  performanee  of  the  COFDM  system  and  pre¬ 
sented  the  BER  performanees  of  the  system  under  four  ehannels.  The  effeets  of  AWGN, 
multipath,  and  AWGN  and  multipath  effeets  over  QPSK  constellations  were  examined. 
As  the  main  of  goal  of  this  thesis,  the  ehapter  also  presented  the  performanee  of  the 
COFDM  system  under  mobile  outdoor  ehannel  eonditions.  The  roles  played  by  RM 
eodes,  peak  elipping,  and  Hanning  windowing  on  PAPR  reduetion  were  also  reported. 
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V.  CONCLUSION 


In  this  thesis,  a  COFDM  based  digital  eommunieation  system  with  Reed-Muller 
error  eorreetion  eoding  was  sueeessfully  simulated.  Three  different  types  of  RM  eodes 
and,  for  eomparison,  a  rate  one-half  eonvolutional  eode  were  investigated  through  simu¬ 
lation  studies  using  QPSK  modulation  type. 

In  order  to  test  the  COFDM  system  in  both  indoor  and  outdoor  ehannel  environ¬ 
ments,  we  ereated  a  step-by-step  test  methodology.  Using  four  different  simulation  steps, 
both  indoor  and  outdoor  ehannel  eharaeteristies  were  studied.  The  bit-error  performanee 
eurves  for  eaeh  of  the  RM  eodes  and  the  eonvolutional  eode  were  obtained. 

Also,  Hanning  windowing  and  peak  elipping  teehniques  were  investigated  for  re- 
dueing  the  PAPR  in  OFDM  systems. 


A.  CONCLUSIONS 

The  results  showed  that  the  COFDM  system  is  robust  in  indoor  ehannel  environ¬ 
ments.  However,  for  the  outdoor  ehannel  environments,  the  system  requires  higher  sig- 
nal-to-noise  ratios  (SNR)  to  aehieve  the  same  bit-error  rate  (BER)  performanee  as  in  the 
ease  of  the  indoor  channel.  Also,  BER  performance  curves  showed  that  the  RM  codes 
provide  better  performance  in  indoor  channel  environment  at  low  SNR.  However,  convo¬ 
lutional  codes  provide  better  performance  in  the  outdoor  channel  environments. 

Reed-Muller  codes  are  straightforward  to  implement,  and  they  provide  a  wide 
range  of  coding  options.  However,  RM  codes  did  not  provide  the  expected  improvement 
in  PAPR  reduction  as  claimed  by  [15].  The  results  showed  that  they  have  almost  the  same 
performance  in  PAPR  reduction  as  convolutional  codes. 

The  addition  of  Hanning  windowing  and  peak  clipping  improve  PAPR  reduction, 
improvement  with  Hanning  windowing  is  limited.  The  results  showed  that,  as  the  width 
of  the  window  is  increased,  the  bit-errors  increased  as  well.  Peak  clipping  provides  con¬ 
siderable  reduction  in  PAPR  but  at  the  cost  of  increased  bit  error.  This  presents  a  tradeoff 
between  the  PAPR  reduction  and  the  BER  performance. 
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B,  FUTURE  WORK 

The  performance  of  a  COFDM  based  digital  communication  system  with  three 
different  Reed-Muller  codes  has  been  successfully  examined  in  this  thesis.  However, 
simulating  the  system  with  different  (r,  m)  pairs  would  be  of  interest.  We  recommend  that 
the  RM  codes  be  further  investigated  to  find  optimal  {r,m)  pairs  to  provide  improvement 
in  PAPR  reduction. 

In  the  simulations,  we  used  QPSK  modulation  technique.  In  future  work,  16-  and 
64-QAM  modulation  techniques  may  be  used.  Therefore,  M-ary  QAM  can  be  examined 
in  a  future  work. 

As  described  in  Chapter  II,  COFDM  based  digital  communication  systems  suffer 
from  frequency  offset  and  are  sensitive  to  time  and  frequency  synchronization.  The  ef¬ 
fects  of  frequency  errors  in  a  COFDM  system  are  analyzed  in  [6].  The  work  in  this  thesis 
and  that  in  [6]  can  be  combined  to  analyze  the  effects  of  Reed-Muller  error  correction 
coding  on  frequency  errors. 

It  is  well  known  that  the  performance  of  a  decoder  is  significantly  increased  with 
soft  decision  decoding.  Therefore,  in  addition  to  the  hard  decision  decoding  technique 
used  in  this  thesis,  a  COFDM  system  can  be  analyzed  with  a  decoder  utilizing  a  soft  deci¬ 
sion  decoding  technique. 
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APPENDIX  A.  SIMULATION-ARCHITECTURE 


The  MATLAB  m-files  used  in  the  simulation  are  deseribed  in  this  ehapter.  In  or¬ 
der  to  aehieve  the  purpose  of  this  thesis,  the  COFDM  simulation  eode  in  [4]  was  modi¬ 
fied  and  adopted.  Sinee  some  of  the  subroutines  in  [4]  were  modified  and  some  were  di- 
reetly  used  in  the  MATLAB  programming,  eitation  to  [4]  was  intentionally  omitted  in  the 
text  in  this  appendix. 


A,  GENERAL  STRUCTURE  OF  MATLAB  PROGRAMMING 

Sinee  the  main  goal  of  this  thesis  was  to  simulate  the  COFDM  system  by  utilizing 
RM  eodes,  all  signal  proeessing  ineluding  forward  error  eorreetion  eoding  (FEC)  and  dif¬ 
ferent  types  of  ehannel  oharaeteristies  were  performed  at  the  baseband. 


1.  COFDM  Transmitter 

The  transmitter  part  of  the  COFDM  system  is  shown  in  Figure  26.  The  subbloeks 
of  the  transmitter  are  explained  in  the  following  diseussion. 


Figure  26.  COFDM  Transmitter  Bloek  Diagram. 


a.  Message  Generator 

This  bloek  generates  a  random  bit  sequenee  representing  the  data  bits.  The 
length  of  the  message  is  variable  and  seleeted  by  the  user.  In  order  to  use  the  same  ran¬ 
dom  message  in  different  simulations,  a  seed  parameter  is  used  in  this  subbloek.  As  long 
as  the  seed  value  remains  the  same,  the  message  generator  generates  the  same  random 
message  and  makes  it  possible  to  run  the  simulation  with  different  parameters  and  eom- 
pare  the  results  for  the  same  information  bit  sequenee. 
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b.  FEC  Coding 

As  described  in  Chapter  III,  in  the  IEEE  802.1  la  standard,  data  are  en¬ 
coded  with  a  convolutional  encoder  with  the  coding  rate  r  =  1/2, 2/3 ,  or  3/4  .  The  convo¬ 
lutional  encoder  uses  the  half-rate  industry-standard  generator  polynomials,  gg  =  133g 
and  gj  =  171g  [7].  In  addition  to  convolutional  codes  used  in  the  simulation  in  [4],  RM 

codes  are  also  used  in  this  thesis.  The  simulation  allows  the  user  to  choose  either  one  of 
these  two  EEC  code  types  or  no  EEC  code.  If  the  convolutional  codes  are  chosen,  the  n 
and  A  of  the  code  are  chosen  by  the  user.  If  the  RM  codes  are  desired,  the  r  and  m  pair  is 
chosen  by  the  user. 


c.  Interleaver 

After  determining  the  length  of  the  message  sequence,  the  simulation  pro¬ 
gram  presents  all  possible  interleaver  pairs  on  the  screen.  According  to  the  interleaver 
dimension  selected  by  the  user,  the  block  interleaver  takes  a  block  of  encoded  bits  and 
permutes  them  as  described  in  Chapter  II. 

d.  Symbol  Reformatter 

Since  binary  phase  shift  keying  (BPSK)  or  quaternary  phase  shift  keying 
(4-PSK  or  QPSK)  can  be  used  in  the  simulation,  the  symbol  reformatter  block  resizes  the 
symbol  lengths  as  1-bit  or  2-bit  in  order  to  form  the  BPSK  or  QPSK  signal  format. 


e.  Differential  Encoder 

Before  the  symbols  are  mapped  into  the  signal,  differential  encoding  is 
performed.  Differential  encoding  is  based  on  a  cumulative  summation.  Cumulative  sum¬ 
mations  are  performed  by  adding  in  the  modulo-A  format.  If  the  frequency  differential 
encoding  is  performed,  then  the  differential  encoding  is  accomplished  along  the  row  vec¬ 
tors  of  the  complex-valued  modulation  matrix.  The  first  element  is  added  to  the  second 
element  in  the  row,  and  the  result  is  replaced  as  the  second  element.  The  first,  second  and 
third  elements  are  added  and  replaced  as  the  third  element  and  this  operation  is  repeated 
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until  all  elements  in  the  row  are  eompleted.  If  the  time  differential  eneoding  is  performed, 
this  operation  is  aeeomplished  along  the  eolumn  veetors  of  the  eomplex-valued  modula¬ 
tion  matrix. 


f.  IFFT  Processing 

The  A^-point  IFFT  operation  is  performed  in  order  to  eonvert  the  frequeney 
domain  data  into  the  time  domain.  Although  A  =  64  is  used  for  all  the  tests  in  this  thesis, 
this  value  ean  be  ehanged  in  the  cofdmsim.m  maero  file. 


g.  Guard  Interval  Insertion 

In  order  to  eliminate  the  ISI  and  ICI  effeets,  16  time  domain  samples  are 
added  to  the  output  of  the  IFFT  proeessor  in  the  IEEE  802.1  la  standard.  In  this  simula¬ 
tion,  this  value  remains  fixed  at  16  as  in  the  standard.  However,  if  it  is  neeessary  to  ex¬ 
tend  this  value  in  order  to  eliminate  the  more  severe  ISI  or  ICI  effects,  this  value  can  be 
changed  in  the  program. 


2.  COFDM  Receiver 

The  operations  in  the  receiver  part  of  the  system  are  the  inverse  operations  of  the 
transmitter.  The  receiver  part  of  the  system  is  shown  in  Eigure  27.  The  subblocks  of  the 
receiver  are  explained  below. 


Received 

baseband 

signal 


Eigure  27.  COEDM  Receiver  Block  Diagram. 


a.  Guard  Interval  Removal 

The  time  domain  samples  added  in  the  transmitter  to  reduce  or  remove  the 
ISI  and  ICI  effects  are  removed  in  this  block  in  order  to  pass  the  actual  information  block 
to  the  sunsequent  subblocks  for  processing. 
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b.  FFT  Processing 

The  A^-point  FFT  operation  is  performed  to  eonvert  the  time  domain  data 
into  the  frequeney  domain.  Due  to  the  multipath  effeets  of  the  transmission  ehannel,  the 
orthogonality  of  the  earners  may  be  distorted.  Therefore,  the  output  of  the  FFT  proeessor 
ean  be  an  array  of  eomplex  values. 

c.  Differential  Decoding  and  Symbol  Reformatting 

Aeeording  to  the  seleetion  of  frequeney  domain  or  time  domain  eneoding, 
the  eorreponding  differential  deeoding  is  performed  in  the  reeeiver.  After  deeoding,  the 
rows  or  eolumns  of  ones  added  by  the  differential  eoding  operation  in  the  transmitter  are 
removed  and  symbols  are  reeonstrueted  aeeording  to  the  magnitude  and  the  phase  of  re- 
eeived  eomplex  values,  whieh  are  the  output  of  the  FFT  proeessor. 

d.  Deinterleaver 

The  output  of  the  symbol  reformatter  is  deinterleaved  to  obtain  the  aetual 
message  order.  Thus,  the  burst  errors  eaused  by  the  noisy  ehannel  are  distributed,  and 
these  randomly  distributed  errors  are  sent  to  the  FEC  deeoder. 


e.  FEC  Decoder 

Depending  on  the  eoding  type  used  in  the  transmitter,  either  the  Viterbi 
deeoder  or  RM  deeoder  deeodes  the  signal  reeeived  from  the  deinterleaver. 

f  Decoded  Message 

The  output  of  the  FEC  deeoder  is  the  reeeived  information  message.  At 
this  point,  depending  on  the  ehannel  eharaeteristies  and  the  eapability  of  the  EEC  deeod¬ 
ing,  some  symbol  errors  may  still  oeeur.  In  order  to  present  the  performanee  of  the  entire 
system,  the  BER  is  ealeulated  and  illustrated  by  the  BER  versus  graphs.  Addi¬ 

tionally,  the  PAPR  values  are  also  ealeulated  and  tabulated. 
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3,  Channel  Models 

In  this  simulation,  five  different  ehannels  are  used.  Channel  0  is  a  noise-free 
ehannel  and  is  performed  to  verify  whether  or  not  the  MATLAB  eoding  works  appropri¬ 
ately.  Channel  1  is  an  AWGN  ehannel.  Channel  2  is  the  multipath  ehannel  that  is  eharae- 
terized  by  frequeney  seleetive  fading  in  dB  and  multipath  delays  in  mieroseeonds.  For 
this  ehannel,  Doppler  shifts  of  5  Hz,  10  Hz  and  15  Hz  were  used.  Additionally,  a  eustom 
ehannel  ean  be  eonfigured  in  aeeordanee  with  the  loss,  frequeney  shifting  and  delay  val¬ 
ues  entered  by  the  user.  Channel  3  is  the  eombination  of  Channels  1  and  2,  i.e.,  it  has 
both  AWGN  and  multipath  ehannel  eharaeteristies.  Channel  4  represents  a  mobile  out¬ 
door  ehannel  that  ineludes  AWGN  and  severe  multipath  effects.  For  this  channel,  two 
sets  of  parameters  are  used  having  different  delays,  fading,  and  Ricean  factors.  The  char¬ 
acteristics  of  the  channels  are  listed  in  Table  8. 


Channel 

Link  Num¬ 
ber 

Noise 

Multipath  El 

ffects 

Fading 

Doppler 

Delay 

#0 

- 

- 

- 

- 

- 

#1 

X 

- 

- 

- 

- 

#2 

1 

- 

X 

5  Hz 

X 

2 

10  Hz 

3 

15  Hz 

4 

Customizable 

#3 

1 

X 

X 

5  Hz 

X 

2 

10  Hz 

3 

15  Hz 

4 

Customizable 

#4 

A 

X 

X 

B 

X 

X 

Tables.  Channel  Characteristics. 


B,  MATLAB  PROGRAMMING  DETAILS 

As  described  in  the  previous  section,  the  simulation  consists  of  three  main  blocks 
which  are  the  receiver,  channel  and  transmitter.  The  main  subroutine  of  the  simulation  is 
chancdl.m,  which  first  calls  cdrcdlft.m  for  the  transmission  processing  and  then  calls  a 
number  of  subroutines  for  the  different  types  of  channels  and  finally  calls  decdrcdl.m  for 
the  receiver  processing.  The  schematic  arrangement  of  the  simulation  is  shown  in  Figure 
28. 
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Figure  28.  General  Arrangement  of  the  COFDM  System. 


In  the  transmitter  part  of  the  system,  message  encoding  is  performed  by  cdrcdlft.m 
and  conversion  from  frequency  domain  to  time  domain  and  guard  interval  insertion  is 
performed  by  tda.m.  The  inverse  operations  are  performed  by  itda.m  and  decdrcdl.m  in 
the  receiver.  The  m-file  itda.m  is  responsible  for  the  guard  interval  removal  and  conver¬ 
sion  to  the  frequency  domain.  The  message  decoding  is  performed  by  decdrcdl.m.  The 
details  of  these  m-files  are  given  in  the  next  subsection. 


1,  COFDM  Transmitter 

The  arrangement  of  the  m-files  that  interact  with  cdrcdlft.m  are  shown  in  Figure 
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Figure  29.  The  Arrangement  of  the  Maero  Files  within  the  edredlft.m. 

A  random  message  of  length  k  is  generated  by  msg.m.  The  length  k  depends  on 
the  number  of  the  symbols  ehosen  by  the  user  and  the  parameter  q,  whieh  represents  the 
number  of  bits  that  eonstruet  eaeh  M-ary  symbol,  where  M  =  2^ .  The  seed  parameter  s 
makes  it  possible  to  test  the  system  with  a  different  parameter  but  the  same  random  mes¬ 
sage  for  eomparing  the  results.  After  the  random  message  sequenee  is  ereated,  it  is  en- 
eoded  by  cnv_encd.m  or  rmencoder.m  depending  on  whether  eonvolutional  eodes  or  RM 
eodes  are  ehosen  by  the  user.  If  the  FEC  eoding  is  not  performed,  cnv_encd.m  is  per¬ 
formed  with  the  same  n  and  k  values.  (Note  that  the  parameter  k  here  is  not  the  length  of 
the  message,  but  represents  the  input  bloek  length  of  the  eonvolutional  eneoder.)  If  RM 
eodes  are  seleeted,  the  monomial  veetors  of  the  RM  eneoding  matrix  are  ereated  by 
gm _partl.m  and  the  remaining  rows  of  the  eneoding  matrix  are  ereated  by  rmgenmat.m, 
using  the  monomial  veetors.  After  that,  the  message  sequenee  is  eoded  by  rmencoder.m. 
If  eonvolutional  eodes  are  seleeted,  the  message  sequenee  is  eoded  by  conv_encd.m.  Af¬ 
ter  FEC  eoding,  the  eoded  bit  sequenee  is  returned  baek  to  marymsg.m,  whieh  then  ealls 
bm.m  to  form  symbols  by  grouping  ^-bits  together  to  ereate  M-ary  symbols.  The  output  is 


49 


a  vector  of  equivalent  decimal  numbers.  The  macro  marymsg.m  calculates  these  values 
using  the  other  subroutines  according  to  the  input  parameters  q,  n,  m,  ecc,  r_rm,  and 
m_rm  where  n  is  the  number  of  the  rows  and  m  is  the  number  of  the  columns  of  the  over¬ 
all  output  message  matrix.  The  parameter  ecc  represents  the  FEC  coding  option,  that  is, 
the  convolutional  codes  or  RM  codes,  and  r_rm  and  m_rm  parameters  are  the  r  and  m 
values  of  the  R{r,  m)  code.  (The  parameter  m  also  represents  the  number  of  the  COFDM 
frequency  tones.  It  must  be  an  even  positive  integer  in  order  to  completely  fill  the  as¬ 
signed  bandwidth.)  After  the  coded  symbol  values  are  returned  back  to  cdrcdlft.m,  this 
function  calls  the  subroutine  cdlilv.m  to  interleave  the  array.  The  subroutine  parameters 
are  I,  k,  dcase,  s  and  sync,  where  /  represents  the  number  of  the  rows  and  k  represents  the 
number  of  the  columns  in  the  interleaver  matrix.  The  parameter  dcase  is  a  variable  indi¬ 
cating  the  interleaving  method  to  be  used.  Although  there  are  nine  cases  in  the  simulation 
program,  only  dcase  0,  which  indicates  the  simple  interleaving,  is  performed  in  this  thesis 
and  the  others  are  beyond  the  scope  of  this  thesis.  The  parameter  s  represents  the  input 
message  to  be  interleaved  and  sync  represents  the  frame  synchronization  bits.  However, 
this  parameter  was  also  not  used  in  this  thesis. 

The  interleaved  message  array  is  converted  from  the  M-ary  to  N-ary  format  for 
the  N-PSK  modulation.  Although  QPSK  is  used  in  this  thesis,  this  value  remains  as  a 
variable  for  different  PSK  modulation  types.  This  conversion  is  conducted  by  bm.m  and 
mb.m.  The  signal  in  the  M-ary  format  is  first  converted  to  binary  by  mb.m  and  then  con¬ 
verted  to  the  N-ary  format  by  bm.m. 

The  message  array  that  is  interleaved  and  converted  from  M-ary  to  N-ary  is  then 
differentially  encoded  since  the  message  array  symbols  are  represented  in  decimal  nota¬ 
tion.  The  transition  of  symbols  to  complex  values  is  performed  by  difcdrft.m.  The  output 
of  this  subroutine  is  an  array  that  represents  the  differentially  encoded  complex  values. 
The  symbol-to-complex  modulation  value  mapping  is  performed  in  two  ways  depending 
on  whether  the  differential  encoding  is  processed  in  time  or  frequency. 

The  last  subroutine  used  within  cdrcdlft.m  is  cmv2fa.m,  which  takes  the  complex 
valued  array  M  and  the  number  of  the  EFT  points  N  as  inputs  and  arranges  M  into  a  fre¬ 
quency  array.  The  output  is  the  shifted  frequency  array,  which  has  padded  zeros  in  the 
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middle  of  the  array.  The  number  of  these  zeros  is  equal  to  the  differenee  between  the 
number  of  FFT  points  and  the  number  of  the  eomplex  values.  Although  RF  eonversion 
and  fdtering  is  not  ineluded  in  this  simulation,  the  padded  zeros  are  eonsidered  as  a  guard 
band  for  the  filter  slopes.  At  this  point,  cdrcdlft.m  eompletes  its  task  and  sends  the  fre- 
queney  array  to  chancdl.m. 

In  order  to  generate  the  COFDM  frequeneies,  the  maero  file  tda.m  is  ealled  by 
chancdl.m.  The  IFFT  operation  is  performed  on  eomplex  modulation  frequeney  array  by 
tda.m,  and  the  periodie  guard  interval  is  also  inserted.  The  input  parameter  is  the 

number  of  the  additional  time  domain  samples  added  to  the  beginning  of  the  information 
symbol  interval  as  deseribed  in  Chapter  II.  The  output  of  this  subroutine  is  the  time  do¬ 
main  samples  prepared  to  be  transmitted  through  the  ehannel. 

The  input  and  output  parameters  of  the  subroutines  that  performed  in  the  transmit¬ 


ter  part  of  the  COFDM  simulation  are  listed  in  Table  9. 


The  name  of 
the  subroutine 

Which  pro¬ 
gram  calls 

Input  parameters 

Output  parameters 

edredlft.m 

chancdl.m 

picy_n,  pic,  dcase,  s, 
freqno,  rintlv,  N, 
mary,nary,fort 

Fa,  MD,  B  ee, 

B  random,  nsymno 

marymsg.m 

cdrcdlft.m 

q,  n,  m,  ecc,  r_rm,  m_rm 

vmary  ce,  random  bit, 
encoutput 

bm.m 

cdrcdlft.m, 

marymsg.m 

q,  V 

m 

msg.m 

marymsg.m 

k 

u 

env  ened.m 

marymsg.m 

ce  g,  ce  kO,ce  input 

ce  output 

rmeneoder.m 

marymsg.m 

r,  m,  msgl 

encd  msg 

rmgenmat.m 

rmeneoder.m 

r,  m 

GM,k 

gm  partl.m 

rmgenmat.m 

m 

bb 

cdlilv.m 

cdrcdlft.m 

1,  k,  dcase,  s,  SYNC 

si 

rotm.m 

cdlilv.m 

V,  m 

vp,  vn 

mb.m 

cdrcdlft.m 

q,  m 

b 

diffedrft.m 

cdrcdlft.m 

q,  m,  fort 

MD 

omv2fa 

cdrcdlft.m 

N,M 

X 

Table  9.  The  Input  and  Output  Parameters  of  the  Subroutines  Run  within  the  Transmitter 

Part  of  the  COFDM  Simulation. 
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2,  COFDM  Receiver 

The  various  m-files  used  in  the  simulations  of  the  COFDM  reeeiver  and  their  in¬ 
teractions  are  shown  in  Figure  30.  The  receiver  operations  are  the  inverse  of  those  in  the 
transmitter. 


Figure  30.  The  Arrangement  of  the  Macro  Files  within  the  decdrcl.m. 


The  complex  modulation  array  is  obtained  from  the  frequency  array  by  using 
fa2cma.m.  The  input  parameters  are  X and  K  where  X is  the  complex  frequency  array  to 
be  rearranged  and  K  is  half  the  number  of  COFDM  frequency  tones.  The  output  parame¬ 
ter  Mn  is  the  complex  modulation  array  of  actual  interleaved  frequencies.  After  the  fre¬ 
quency  array  is  reconstructed  into  the  complex  array,  the  complex  modulation  values  are 
differentially  decoded  into  N-ary  symbols  by  dfdcdrft.m.  All  reference  ones  added  in  the 
transmitter  are  removed,  and  the  remaining  received  message  matrix  is  decoded.  The  in¬ 
put  parameters  MD,  fort,  q  and  qp,  where  MD  is  the  complex  modulation  values  as  the 
output  of  the  m-file  fa2cma.m  and  q  is  the  number  of  bits  in  a  symbol  according  to  the  N- 
ary  format.  The  parameter  fort  represents  whether  the  differential  coding  is  performed  in 
frequency  or  time.  The  output  parameters  of  this  m-file  are  s  and  M,  where  s  represents 
the  corresponding  mapped  symbols  in  decimal  format  and  M  is  the  differentially  decoded 
modulation  array.  In  order  to  convert  the  modulation  array  values  from  N-ary  to  M-ary, 
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values  are  eonverted  from  N-ary  to  binary  by  mb.m  and  binary  to  M-ary  by  bm.m.  The  m- 
file  cdldlv  deinterleaves.  The  input  parameter  si  represents  the  reeeived  interleaved  mes¬ 
sage  array  and  the  other  parameters  are  the  same  on  the  ones  in  cdlilv.m.  The  m-file 
rotm.m  is  ealled  by  cdlilv.m  to  rotate  the  array. 

The  input  and  output  parameters  of  the  subroutines  used  in  the  reeeiver  part  of  the 


COFDM  simulation  are  listed  in  Table  10. 


The  name  of  the 
subroutine 

Which  program 
calls 

Input  parameters 

Output  parameters 

decdrcdl.m 

Chancdl.m 

picy  n,  pic,  dcase,  K,  Fa, 
nsymno,  ffeqno,  rdintlv, 
cdintlv,  mary,  nary,  fort, 

B  random,  ecc,  r_rm, 
m  rm 

outmsg,  de- 
coder_output,_bit,  ran- 
dom  msg,  random  bit, 

M,  MM,  binary  value 

fa2cma.m 

decdrcdl.m 

K,X 

Mm 

dfdcdrft.m 

decdrcdl.m 

qp,  q,  MD,  fort 

s,  M 

bm.m 

decdrcdl.m 

q,  V 

M 

cdldlv.m 

decdrcdl.m 

1,  k,  dcase,  si,  SYNC 

s 

rotm.m 

cdldlv.m 

V,  m 

vp,  vn 

mb.m 

decdrcdl.m 

q,  m 

b 

viterbi.m 

decdrcdl.m 

v_G,  v_k,  channel  output 

decoder  output,  survi- 
vor  state,  cumu¬ 
lated  metric 

nxt  stat.m 

viterbi.m 

current  state,  input,  v  L, 

V  k 

next  state,  mem¬ 
ory  contents 

bin2deci.m 

nxt  stat.m, 
viterbi.m 

v_x 

v_y 

deci2bin.m 

nxt  stat.m, 
viterbi.m 

X,  1 

y 

rmdecoder.m 

decdrcdl.m 

Me,  r,  m 

Msg  decoded  line 

rmgenmat.m 

rmdecoder.m 

r,  m 

GM,k 

gm  parti 

rmgenmat.m 

m 

bb 

WW  mat2.m 

rmdecoder.m 

GM,  r,  m  ,k 

WWl 

inverse.m 

WW  mat2.m 

D 

d 

dot_pro.m 

rmdecoder.m 

X,  S 

z 

Table  10.  The  Input  and  Output  Parameters  of  the  Subroutines  Run  within  the  Reeeiver  Part 

of  the  COFDM  Simulation. 


3.  COFDM  Channel 

Figure  3 1  shows  the  m-files  used  to  simulate  the  ehannel. 
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TRANSMITTER 


cdrcdift.m 


tda.m  * 


CHANNEL  MODEL  1 


awgn.m  ♦ 


chancdl.m 


CHANNEL  MODEL  2 


*|  chuhf.m 


raydop.m 


ofst.m 


H  dline.m 

-f— 


cvdd.m 


CHANNEL  MODEL  3 


check. m 


RECEIVER 


decrcdl.m 


►  itda.m 


’ 

f 

channeH  .m 

jakes.m 

CHANNEL  MODEL  4 

itda.m 


Figure  3 1 .  The  Arrangement  of  the  Maero  Files  within  the  ehanedl.m. 


The  input  parameters  of  the  macro  file  chuhf.m  are  5,  x,  loss,  dly,  dop,  N  and 
freqspace,  where  loss  is  the  signal  loss  in  dB,  dly  is  the  time  delays  in  microseconds  and 
dop  is  the  Doppler  frequency  shifting  in  Hz.  Multipath  delay  effects  are  created  by 
dline.m  and  cvdd.m.  The  input  array  is  filtered  by  a  FIR  filter  whose  coefficients  are  cre¬ 
ated  by  cvdd.m  according  to  the  values  in  the  array  dly.  As  a  next  step,  the  Doppler  shift 
frequency  is  calculated  by  ray_dop.m  according  to  the  seed  value  s,  freqspace  and  the 
number  of  the  FFT  points  N.  The  channel  offset  is  calculated  by  ofst.m.  Finally,  the 
power  losses  are  calculated  by  multiplying  each  loss  value  in  the  array  loss  with  the  cor¬ 
responding  message  value.  The  output  of  chuhf.m  is  a  complex  valued  time  domain  array 
including  multipath  effects. 

Channel  3  is  realized  by  using  awgn.m  and  chuhf.m.  Channel  4  is  simulated  by 
channell.m  and  Jakes. m.  The  first  subroutine  calls  the  second  and  provides  two  severe 
channels  having  different  multipath  effects.  The  selection  of  these  two  channels  must  be 
conducted  before  running  the  simulation. 

The  input  and  output  parameters  of  the  subroutines  simulating  the  channel  part  of 
the  COFDM  simulation  are  shown  in  Table  11. 
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The  name  of 
the  subroutine 

Which  pro¬ 
gram  calls 

Input  parameters 

Output  pa¬ 
rameters 

awgn.m 

ehanedl.m 

X,  N,  sigma 

Y 

ehuhfm 

ehanedl.m 

s,  X,  loss,  dly,  dop,  N, 
freqspaee 

y 

raydop.m 

ehuhfm 

s,  M,  N,  es 

e 

ofst.m 

ehuhfm 

e,  N,  X 

xo 

dline.m 

ehuhfm 

X,  d 

xd 

evdd.m 

dline.m 

X,  alpha 

y 

ehannell.m 

ehanedl.m 

xt 

yr 

jakes.m 

ehannell.m 

Fd,  Fs,  N 

8 

itda.m 

ehanedl.m 

Ng,y 

y 

Table  1 1 .  The  Input  and  Output  Parameters  of  the  Subroutines  Run  within  the  Channel  Part 

of  the  COFDM  Simulation. 
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APPENDIX  B.  COFDM  MATLAB  CODE 


This  appendix  presents  the  MATLAB  m-files  used  in  the  simulations.  Deserip- 
tions  of  the  m-files  here  ean  be  found  in  Appendix  A. 


AWGN 

%AWGN 

%TITLE  :  ADDITIVE  WHITE  GAUSSIAN  NOISE 

% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  :  TEZEREN,Serdar  Umit,ETJG, TURKISH  NAVY 

% - 

funetion  Y  =  awgn(X,N, sigma) 

%Pind  dimensions  of  the  input  array 
[rr,ee]=size(X); 

%Generate  a  random  real  part 
seed=14; 

randn('state',seed); 

wreal=randn(rr,ee); 

%Generate  a  random  imaginary  part 
seed=seed+6; 
randn('state',seed); 
wimg=i*randn(rr,oo); 

%An  array  of  random  eomplex  entries  ehosen  from  a  normal  distribution  with 
%mean  0.0  and  varianee  1.0.  Array  dimension  is  the  same  as  X. 
W=wreal+wimg; 

%Random  noise  multiplied  by  the  sigma  faetor  and  added  to  the  signal. 
Y=X+(sigma.  *  W); 

% - 
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BIN2DECI 


%BIN2DECI 

%TITLE  :  BINARY  TO  DECIMAL  CONVERSION 

% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  v_y=bin2deoi(v_x) 

v_l=length(v_x); 

v_y=(v_l-l:-l:0); 

v_y=2Av_y; 

v_y=v_x*v_y'; 

% - 
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BM 

%BM 

%TITLE  :  BINARY  TO  M-ARY  CONVERTER 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  m=bm(q,v) 

%Pind  the  length  of  input  veotor,v,and  determine  if  there  is  a  remainder 

%after  dividing  by  q 

n=length(v); 

r=rem(n,q); 

%If  there  is  no  remainder, don't  pad  v  input  veetor.  Otherwise  add  the  appropriate 

%number  of  zeros  to  generate  a  eode  word  with  an  exaet  multiple  of  q  bits. 

ifr==0 

v=v; 

else 

v=[v  zeros(l,q-r)]; 
end 

%Plaee  least  signilieant  bit  of  the  symbol  on  the  left  end. 

map=l; 

for  j=l;q-l 

map=[map  2^j]; 

end 

%Remove  q  bits  at  a  time  from  v  to  generate  m-ary  symbol  values. 

n=length(v); 

p=round(n/q); 

A=zeros(q,p); 

A(:)=v; 

m=map*A; 

m_ary_msg=m; 

% - 
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CDLDLV 


%CDLDLV 

%TITLE  :  BLOCK  DEINTERLEAVER 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof.  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  s=cdldlv(l,k,doase,si,SYNC) 

si(length(si)+l-length(SYNC):length(si))=zeros(l,length(SYNC)); 

N=length(si); 

ifl*k==N 

x=zeros(l,k); 

x(:)=si; 

KKl:k)-l; 

CR=K.*(K+l)/2; 

RR=L.*(L+l)/2; 
if  doase==l 
for  kk=l;k 

x( :  ,kk)=rotm(x( :  ,kk),  CR(kk)) ; 
end 

elseif  dease==2 
for  kk=l;k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

elseif  dease==3 
for  kk=l;l 

x(kk,  :)=rotm(x(kk,  :),RR(kk)); 
end 

elseif  dease==4 
for  kk=l;l 

[z,x(kk,:)]=rotm(x(kk,:),RR(kk)); 

end 

elseif  dease==5 
for  kk=l:k 

x( :  ,kk)=rotm(x( :  ,kk),  CR(kk)) ; 
end 

for  11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

elseif  dease==6 
for  kk=l:k 
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[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

for  11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

elseif  dcase==7 
for  kk=l:k 

x( :  ,kk)=rotm(x( :  ,kk),  CR(kk)) ; 
end 

for  11=1:1 

[z,x(ll,:)]=rotm(x(ll,:),RR(ll)); 

end 

elseif  dcase==8 
for  kk=l:k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

for  11=1:1 

[z,x(ll,:)]=rotm(x(ll,:),RR(ll)); 

end 

end 

x=x'; 

s=x(:); 

s=s'; 

end 

% - 
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CDLILV 


%CDLILV 

%Title  :  BLOCK  INTERLEAVER 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof.  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  si  =  odlilv(l,k,doase,s,SYNC) 

N=length(s); 

ifl*k==N 

x=zeros(l,k); 

x=x'; 

x(:)=s; 

x=x'; 

Intermediate_mx=x; 

K=(l:k)-1; 

CR=K.*(K+l)/2; 

RR=E.*(E+l)/2; 
if  doase==l 
for  kk=l:k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

elseif  dease==2 
for  kk=l;k 

x( :  ,kk)=rotm(x( :  ,kk),  CR(kk)) ; 
end 

elseif  dease==3 
for  kk=l;l 

[z,x(kk,:)]=rotm(x(kk,:),RR(kk)); 

end 

elseif  dease==4 
for  kk=l;l 

x(kk,  :)=rotm(x(kk,  :),RR(kk)); 
end 

elseif  dease==5 
for  11=1:1 

[z,x(ll,:)]=rotm(x(ll,:),RR(ll)); 

end 

for  kk=l:k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 
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elseif  dcase==6 
for  11=1:1 

[z,x(ll,:)]=rotm(x(ll,:),RR(ll)); 

end 

for  kk=l:k 

x( :  ,kk)=rotm(x( :  ,kk),  CR(kk)) ; 
end 

elseif  dease==7 
for  11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

for  kk=l:k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

elseif  dease==8 
for  11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

for  kk=l:k 

x( :  ,kk)=rotm(x( :  ,kk),  CR(kk)) ; 

end 

end 

si=x(:); 

si=si'; 

end 

si(length(si)-length(SYNC)+ 1  :length(si))=SYNC; 
% - 
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CDRCDLFT 


%CDRCDLFT 

%Title  :  COFDM  ENCODER  WITH  BEOCK INTEREEAVER 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof.  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  :  TEZEREN,Serdar  Umit,ETJG, TURKISH  NAVY 

% - 

fune- 

tion[Pa,MD,B_ee,B_random,nsymno,ene_output]=edredlft(piey_n,pie,dease,s,freqno,rint 
lv,eintlv,N,mary, nary, fort, eee,r_rm,m_rm); 

%  Determine  whether  the  number  of  OEDM  frequeneies  are  even  (#  of  matrixeol- 
umns),indioated 

%  by  the  "fireqno"  parameter.  If  odd  go  to  error  message.  Odd  frequeneies  are  not  allowed 
%  sinee  the  formation  of  the  frequeney  array  is  symmetrieal  and  even, 
if  rem(freqno,2)~=0 

disp('ERROR;  The  number  of  matrix  oolumns,freqno,representing  OEDM  frequeneies 
must  be  an  even  number!') 
elseif  rem(freqno,2)==0 

%  Determine  if  the  row  and  eolumn  interleave  parameters  are  greater  than  fireqno  when 
%  multiplied  together.  If  not,  then  display  error  message  and  stop, 
if  (rintlv*ointlv)<(freqno) 
disp(") 

disp('ERROR;  The  row  and  eolumn  interleave  parameters  are  not  eompatible  with  #  of 

OEDM  frequeneies!') 

dispC) 

else 

%  Caleulate  the  row  symbol  number 
symno=rintlv  *  eintlv/ fireqno ; 

%  Display  error  message  if  symno  and  fireqno  not  eompatible  with  rintlv  and  eintlv  and 
stop. 

%  If  not  oompatible,the  interleaver  funetion  does  not  work  eorreetly. 

if  rem(symno,l)~=0 

dispC) 

disp('ERROR:  The  row  and  eolumn  interleave  parameters  are  not  eompatible  with  #  of 
OEDM  frequeneies!') 

dispC  Eor  the  entered  rintlv,  eintlv,  and  fireqno  parameters,  the  ealeulated  symno  is:') 
disp(symno) 

multiesall=mltpl(rintlv,ointlv); 

multies=multiesall(  1 ,  (2 :  length(multiesall)- 1 )) ; 

dispC  Possible  ehoiees  for  fireqno  based  upon  rintlv  and  eintlv  are:') 

dispC) 

disp(multies) 
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elseif  rem(symno,l)==0 
if  freqno  >=  N; 
disp(") 

disp('ERROR;  The  number  of  frequeney  points,  N,  needs  to  be  inereased  !') 

disp('N  must  be  larger  than;') 

disp(") 

disp(freqno) 

disp(") 

elseif  freqno  <  N; 

%  Generate  a  random  message  matrix  of  m-ary  symbols, based  upon  parameter,mary. 
Nmbrofsymbols  =symno*  freqno; 

[B_oe,B_random,eno_output]=marymsg(mary,symno,freqno,eoo,r_rm,m_rm); 

Rndm_m_ary_msg=B_random; 

%  Perform  a  bloek  interleaving  funetion  on  the  matrix,  B,  with  rintlv  rows 
%  and  eintlv  eolumns. 

SYNC=[]; 

[Br  Bo]=size(B_oe); 

Bt=B_oe'; 

Bveot=Bt(:)'; 

si=odlilv(rintlv,ointlv,doase,Bveot,SYNC); 

Bi=reshape(si,Bo,Br)'; 

Intrlvd_array=Bi; 
ml=bm(nary,mb(mary,Bi)); 
lengthm  1 =length(m  1 ) ; 
nsymno=lengthml ; 
remml=rem(lengthml , freqno); 
if  remml==0; 
ml=ml; 
else 

zero=zeros(freqno-remm  1 ); 

ml=[ml  zero(l,:)]; 

end 

length2m  1 =length(m  1 ) ; 

m=(reshape(m  1 ,  freqno ,  length2m  1  /  freqno))' ; 

N_ary_msg=m; 

%  Generate  a  differentially  eneoded  matrix  of  eomplex 
%  values  with  unit  magnitude  and  one  of  (2^n)  equal  phases. 
MDD=difedrft(nary,m,fort); 

[MDm  MDn]=size(MDD); 

MD=MDD ;  Cmplx_mod_array=MDD ; 

%  Form  the  frequeney  array  of  modulation  values  that  inelude  guard  interval. 
Fa=emv2fa(N,MD); 

Freq_array=Fa; 
end;  end 
end;  end 


%■ 
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CHANCDL 


%CHANCDL 

%  TITLE  :  SIMULATIONS  FOR  AWGN  &  MULTIPATH  FADING  CHANNEL 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Dave  Roderiek,  Naval  Postgraduate  Sehool 

%  Modified  by  :  TEZEREN,Serdar  Umit,LTJG, TURKISH  NAVY 

% - 

fune- 

tion[errmax, errors, freqerrs,papr]=ehanedl(ehnmdl, wait, prut, piey_n, pie, dcase,s,freqno,rint 

lv,eintlv,N,mary,nary,n,k,blklgth,Ng,sigs,loss,dly,dop,freqspaee,fort,eee,r_rm,m_rm) 

sigveet=sigs; 

klgth=length(k); 

ehklp=l; 

errveet=[]; 

berveet=[]; 

freqerrmx=[]; 

errsperpr=[]; 

Es_No=[]; 

Eb_No=[]; 

sermx=[]; 

bermx=[]; 

rowerrmx=[]; 

symno=rintlv  *  eintlv/ freqno ; 
for  lp=l;length(sigveot); 

% 

[xmt,modvals,B_oe,B_random,nsymno,eno_output]=odrodlft(pioy_n,pio,doase,s,freqno,ri 

ntlv,ointlv,N,mary,nary,fort,eoo,r_rm,m_rm); 

xmtifft=tda(N  g  ,xmt) ; 

[ppp]=papr3(xmtifft,N,freqno,pio); 

xmtifft=ppp; 

size_modvals=size(modvals) 

xmtpts=  1  ;length(xmtifft); 
if  ohnmdl==0 
sandn=xmtifft; 
elseif  chnmdl==1 

disp(['Sigma=',nuni2str(sigveot(lp))]); 
sandn=awgn(xmtifft,N,sigveot(lp)); 
elseif  ehnmdl==2 

sandn=ehuhf(s+l,xmtifft,loss,dly,dop,N,freqspaee); 
elseif  ehnmdl==3 

sandmltpth=ehuhf(s+l,xmtifft,loss,dly,dop,N,freqspaee); 
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disp(['Sigma=',nuni2str(sigvect(lp))]); 
sandn=awgn(sandmltpth,N,sigvect(lp)); 
elseif  chnmdl==4 
sandnmlt=channel  1  (xmtifft); 
sandn=awgn(sandnmlt,N,sigvect(lp)); 
end 

sandnfft=itda(Ng,sandn); 

K=(length(modvals(l  ,:)))/2; 

[rcvd,rcvd_bit,random_msg,random_bit,M,MM,binary_value]=decdrcdl(picy_n,pic,dcase 
,K,sandnfft,nsymno,freqno,rintlv,cintlv,mary, nary, fort, B_random,ecc,r_rm,m_rm); 
Transmitted_msg=random_msg; 

Received_msg=rcvd; 

[er- 

rors,bit_error,freqerrs,errmx,rowerrs]=check(pic,random_msg,random_bit,rcvd,rcvd_bit, 

n,k(chklp),blklgth,ecc,r_rm,m_rm,enc_output,binary_value); 

errvect=  [errvect,errors] ; 

bervect=[bervect,bit_error]; 

freqerrmx=  [freqerrmx;freqerrs] ; 

rowerrmx=[rowerrmx;rowerrs] ; 

cmtEs_No=  1/(2  *N*  (sigvect(lp)^2)); 

cmtEb_N  o=crntEs_N  o ; 

Es_No=[Es_No,crntEs_No]; 

Eb_N  0= [Eb_N  o  ,crntEb_No] ; 

Es_Nodb=10*logl0(Es_No); 

Eb_Nodb=10*logl0(Eb_No); 

end 

ser=errveet/(symno  *  freqno) ; 

ber=berveet/ (2  *  symno  *  freqno); 

sermx=  [sermx;ser] ; 

bermx=[bermx;ber]; 

errsum=sum(errveot); 

errsperpr=[errsperpr,errsum] ; 

errmax=max(rowerrmx') ; 

if  pioy_n==l 

figure(pio+l) 

plot(modvals,'*') 

hold  on; 

plot(0,0,'+') 

hold  off; 

title(['Transmitted  Signal  ',int2str(2^nary),  '-ary  Constellation  Plot']) 

xlabel(['Magnitude=r]) 

axis('square'); 

orient  tall 
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grid 

if  pmt==l; 
print 

pause(lO); 

end 

pause(wait); 

end 

% 


if  piey_n==l 
figure(pie+2) 
plot([0;N-l],abs(xmt),'*') 

title(['Frequeney  Array  Plot  (The  number  of  FFT  points  are  ',int2str(N),')']) 

xlabel( ['Guard  interval  length  is  ',mt2str(Ng)]) 

axis('square'); 

orient  tall 

grid 

if  pmt==l; 
print 

pause(lO) 

end 

pause(wait); 

end 

% 

if  piey_n==l 
figure(pie+3) 
surf(abs(modvals)); 
shading  interp 
grid 

orient  tall 

title( ['Magnitude  of  Transmitted  Signal']) 
xlabel('OFDM  Freq  #') 
ylabel('Symbol  Row  Number') 
zlabel(  ['Magnitude']) 
if  pmt==l; 
print 

pause(lO) 

end 

%pause(wait); 

end 

% 

if  piey_n==l 
figure(pie+6) 
plot(M,'*') 
hold  on; 
plot(0,0,'+') 
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hold  off; 

title( ['Received  int2str(2^nary),'-ary  Signal  Constellation  Plot,  BEFORE  Differential 

Decoding']) 

orient  tall 

axis('square'); 

grid 

if  pmt==l; 
print 

pause(lO) 

end 

pause(wait); 

end 

% 

if  picy_n==l 
figure(pic+7) 
plot(MM,'+') 
hold  on; 
plot(0,0,'+') 
hold  off; 

title( ['Received  ',int2str(2^nary),'-ary  Signal  Constellation  Plot,  AFTER  Differential  De¬ 
coding']) 
orient  tall 
axis('square'); 
grid 

if  pmt==l; 
print 

pause(lO) 

end 

pause(wait); 

end 

% 

% 

if  picy_n==l 

%roty_n=input('Do  you  want  to  rotate  3-D  plot?(yes=l,no=0):'); 

roty_n=0; 

figure(pic+8) 

surf(abs(M)); 

grid  on 

shading  interp 

%grid 

orient  tall 

title( ['Magnitude  Variation  of  Received  Signal']) 
xlabel('OFDM  Freq  #') 
ylabel('Symbol  Row  Number') 
zlabel(  ['Magnitude']) 
if  roty_n==l 
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%Rotate  the  3-D  plot 
for  k=l:5 

view(-70+10*k,15+10*k) 

dispC); 

disp('Press  "enter"  to  rotate  plot...'); 

pause 

end 

%elseif  roty_n==0 
end 

if  pmt==l; 
print 

pause(lO) 

end 

%pause(wait); 

end 

% 

if  errsum~=0 

%2-D  Error  Performanee  Curve  showing  BER  vs.  Es/No. 

figure(pio+12) 

semilogy(Eb_Nodb,ber) 

grid 

if  fort==l 

if  dop==[5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5] 

title(['Eink  1  :  Performanee  graph;  BIT  Error  Rate  vs.  Eb/No  (Ereq.  Diff  Ene.)(Total 
errors=',int2str(sum(errveet)),')']) 

elseifdop==[10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10] 
title(['Eink  2  :  Performanee  graph:  BIT  Error  Rate  vs.  Eb/No  (Ereq.  Diff.Ene.)(Total 
errors=',mt2str(sum(errveet)),')']) 

elseifdop==[I5,I5,I5,I5,I5,I5,I5,I5,I5,I5,I5,I5,I5,I5,I5,I5,I5,I5] 
title(['Eink  3  :  Performanee  graph;  BIT  Error  Rate  vs.  Eb/No  (Ereq.  Diff.Ene.)(Total 
errors=',mt2str(sum(errveet)),')']) 
else 

title( ['Custom  Eink  Performanee  graph;  BIT  Error  Rate  vs.  Eb/No  (Ereq. 
Diff.Ene.)(Total  errors=',mt2str(sum(errveet)),')']) 
end 

elseif  fort==0 

if  dop==[5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5] 

title(['Emk  I  :  Performanee  graph;  BIT  Error  Rate  vs.  Eb/No  (Ereq.  Diff  Ene.)(Total 
errors=',mt2str(sum(errveet)),')']) 

elseif  dop==[I0, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] 
title(['Eink  2  :  Performanee  graph:  BIT  Error  Rate  vs.  Eb/No  (Ereq.  Diff.Ene.)(Total 
errors=',mt2str(sum(errveet)),')']) 

elseif  dop==[I5, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15] 
title(['Eink  3  :  Performanee  graph:  BIT  Error  Rate  vs.  Eb/No  (Ereq.  Diff  Ene.)(Total 
errors=',mt2str(sum(errveet)),')']) 
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else 

title( ['Custom  Link  Performanee  graph:  BIT  Error  Rate  vs.  Eb/No  (Ereq. 
Diff.Ene.)(Total  errors=',int2str(sum(errveet)),')']) 
end 
end 

xlabel(['Eb/No(dB)(#  of 

OEDM=',int2str(freqno),')(doase=',mt2str(doase),')(Interleaverpair=',int2str(rintlv),',',mt2s 
tr(ointlv),')  M-ary=',int2str(2^mary),',Nary=',int2str(2^nary)]); 

ylabel(['Sigma  Range:(',num2str(min(sigs)),'- 
',num2str(max(sigs)),')(RS=',mt2str(floor((n- 
k)/2)),')(Symbol#=',mt2str(symno*freqno),')(Seed=',num2str(s), ')']); 
orient  tall 

%2-D  Error  Performanee  Curve  showing  SER  vs.  Eb/No. 
figure(pio+13) 
semilogy(Es_Nodb,ser) 
grid 

if  fort==l 

if  dop==[5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5] 

title(['Eink  1  :  Performanee  graph:  Symbol  Error  Rate  vs.  Es/No  (Ereq. 
Diff.Ene.)(Total  errors=',mt2str(sum(errveet)),')']) 

elseifdop==[10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10] 
title(['Eink  2  :  Performanee  graph:  Symbol  Error  Rate  vs.  Es/No  (Ereq. 
Diff.Ene.)(Total  errors=',mt2str(sum(errveet)),')']) 

elseifdop==[15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15] 
title(['Eink  3  :  Performanee  graph:  Symbol  Error  Rate  vs.  Es/No  (Ereq. 
Diff.Ene.)(Total  errors=',mt2str(sum(errveet)),')']) 
else 

title( ['Custom  Eink  Performanee  graph:  Symbol  Error  Rate  vs.  Es/No  (Ereq. 
Diff.Ene.)(Total  errors=',mt2str(sum(errveet)),')']) 
end 

elseif  fort==0 

if  dop==[5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5] 

title(['Eink  1  :  Performanee  graph:  Symbol  Error  Rate  vs.  Es/No  (Ereq. 
Diff.Ene.)(Total  errors=',mt2str(sum(errveet)),')']) 

elseif  dop==[10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] 
title(['Eink  2  :  Performanee  graph:  Symbol  Error  Rate  vs.  Es/No  (Ereq. 
Diff.Ene.)(Total  errors=',int2str(sum(errveet)),')']) 

elseif  dop==[15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15] 
title(['Eink  3  :  Performance  graph:  Symbol  Error  Rate  vs.  Es/No  (Ereq. 
Diff.Enc.)(Total  errors=',int2str(sum(errvect)),')']) 
else 

title( ['Custom  Eink  Performance  graph:  Symbol  Error  Rate  vs.  Es/No  (Ereq. 
Diff.Enc.)(Total  errors=',int2str(sum(errvect)),')']) 
end 
end 

text(min(ceil(Es_Nodb)),.l  8, ['Eoss=[',num2str(loss), ']']); 
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text(min(ceil(Es_Nodb)),.  12, ['Delay=[',num2str(dly), ']']); 

text(min(ceil(Es_Nodb)),. 08, ['Doppler=[',num2str(dop), ']']); 

xlabel(['Es/No(dB)(#  of  OEDM=',int2str(freqno),')(dcase=',int2str(dcase),')(Interleaver 
pair=',int2str(rintlv),',',int2str(cintlv),')  M-ary=',int2str(2^mary),',Nary=',int2str(2^nary)]); 

ylabel(['Sigma  Range:(',num2str(min(sigs)),'- 
',num2str(max(sigs)),')(RS=',int2str(floor((n- 
k)/2)),')(Symbol#=',int2str(symno*freqno),')(Seed=',num2str(s),')']); 

orient  tall 
end 

if  pmt==l 
print 

pause(lO) 

end 

% - 
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CHANNELl 


%CHANNEL1 

%  TITLE  :  SEVERE  CHANNEL 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof  R.  Cristi,  Naval  Postgraduate  Sehool 

%  Modified  by  :  TEZEREN,Serdar  Umit,LTJG, TURKISH  NAVY 

% - 

funetion  yr=ehanneH(xt) 

[asd  asn]=size(xt); 

xt=xt.'; 

xt=xt(:).’; 

%  yr=ehannel(xt) 

%  with  xt=transmitted  veetor 
%  yr=reeeived  veetor 
%  length(yr)=length(xt) 

%  if  the  program  gives  you  an  error,  just  inerease  the  size  of  the 
%  transmitted  veetor  "xt" 

%  Created  by:  Roberto  Cristi,  November  2003 

Ps=l .25*  10^6;  %  sampling  frequeney  (Hz) 

Pd=200;  %  doppler  frequeney  (Hz) 

%  Parameters  for  Channel  A 

%Td=[0.0,  0.25,0.5,1.0,1.9,2.2];  %  time  delays  (in  mierosee) 

%PdB=[0.0,-1.0,-9.0,-10.0,-15.0,-20.0];  %  Powers  (in  dB) 
%K=[0,0,0,0,0,0];  %  Rieean  Eaetor 

%  Parameters  for  Channel  B 

Td=[0.0,  0.25,9.0,13.0,  17.0,  20.0];  %  time  delays  (in  mierosee) 

PdB=[-2.5,0.0,-12.8,-10,  -25.2,  -16];  %  Powers  (in  dB) 
K=[0.5,0.5,0,0,0,0];  %  Rieean  Eaetors 

Td=Td*(10^(-6));  %  time  delays  (in  see) 

nd=round(Td*Ps);  %  time  delay  in  samples 
P=10.^(PdB/10);  %  Powers  (Linear) 

Np=length(xt); 

yr=zeros(size(xt)); 
for  k=  1 :  length(T  d) 
g=jakes(Ed,  Es,  Np); 

s=sqrt(P(k)/(K(k)-l-l));  %  random  path 

m=sqrt(P(k)*K(k)/(K(k)+l));  %  direet  path 
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total=s*g+m*exp(j*2*pi*Fd/Fs*(0;length(xt)-l)); 
yr=yr+[zeros(  1  ,nd(k)),  xt(  1  ;length(xt)-nd(k))] .  *total; 
end 

yr=reshape(yr,asn,asd).'; 

% - 
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CHECK 


%CHECK 

%  TITLE  :  SOURCE  AND  SINK  MESSAGE  CHECKER 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  :  TEZEREN,Serdar  Umit,LTJG, TURKISH  NAVY 

% - 

fune- 

tion[error_no,bit_error_total,freqerrs,errmx,rowerrs]=eheek(pie,x,xbit,y,ybit,n,k,blklgth,e 

ee,r_rm,m_rm,ene_output,binary_value) 

if  eee~=2;%  unless  reed  muller 

if  blklgth>n 

disp(") 

disp('ERROR!  The  bloek  length,  blklgth,  must  be  equal  or  less  than  the  eode  word 
length,n.') 

disp('Please  enter  a  smaller  value  for  blklgth,  or  ehange  n.') 
disp(") 

elseif  blklgth<=n 
if  n<k 
disp(") 

disp('Error!  The  eode  word  length,n,must  be  equal  or  larger  than  the  information 
length,k.') 

disp('Please  enter  a  larger  value  for  n,  or  ehange  k  to  a  smaller  number.') 

disp(") 

elseif  n>=k 

Pirst_matrix=x; 

Seeond_matrix=y; 

[rx  ex]=size(x); 

%Compare  inputs  x  and  y  and  generate  error  matrix,  "errors" 
errors=(x~=y); 

Pirst=xbit; 

Seoond=ybit; 

[rxl  oxl]=size(xbit); 

%Compare  inputs  xbit  and  ybit  and  generate  BIT  error  matrix,  "bit_errors" 
bit_errors=(xbit~=ybit); 

%Pind  the  error  distribution  vs.  OEDM  frequeneies 
freqerrs=sum(  errors); 

%Pind  the  error  loeation  in  "errors"  where  element  in  x  and  y  differ. 
Error_loeations=(lind(  errors))'; 

Error_number=sum(sum(  errors)); 

Correot_symbl_num=(size(y,l)*size(y,2))-Error_number; 


75 


%Find  the  bit  error  location  in  "errors"  where  element  in  xl  and  yl  differ. 
%bit_Error_locations=(lind(bit_errors))'; 
bit_error_total=sum(sum(bit_errors)); 
symcorr=floor((n-k)/2); 
if  blklgth<=(n-k) 

disp('Error! !  !The  block  length  is  too  short  for  the  given  n  and  k  values') 
disp(") 

elseif  blklgth>(n-k) 
errtrans=errors'; 

%Reshape  the  error  matrix  as  a  vector  of  errors 
errvect=errtrans(:)'; 
b  lkrem=rem(length(errvect)  ,blklgth) ; 
if  blkrem~=0; 

zeropad=zeros(blklgth-blkrem); 
errvectpad=[errvect  zeropad(l,:)]; 
elseif  blkrem==0; 
errvectpad=errvect; 
end 

blknos=length(errvectpad)/blklgth; 

errcorct=[]; 

errblksum=[]; 

for  lp=l:blknos; 

errblk=errvectpad(((blklgth*(lp-l))+l):(blklgth*lp)); 

errblklgth=length(errblk); 

if  sum(errblk)<=symcorr; 

noerr=zeros(errblklgth); 

errblk=noerr(l,:); 

elseif  sum(errblk)>symcorr; 

errblk=errblk; 

end 

errcorct=[errcorct  errblk]; 
errblksum=[errblksum  sum(errblk)]; 
end 

newerrvect=errcorct(  1  :length(errvect)); 
errtot=sum(newerrvect) ; 

RSerrs=(reshape(newerrvect,size(errors,2),size(errors,l)))'; 

%Eind  the  error  distribution  vs.  OEDM  Erequencies 
freqerrs=sum(RS  errs) ; 
errindex=(fmd(RSerrs))'; 

RSerrtot=sum(errblksum); 

RS  errdif=Error_number-RS  errtot; 
errperblk=[(l:blknos);errblksum]; 

%Check  to  see  if  x  and  y  are  the  same.  If  not,  display  error  message 
if  x==y; 

disp('GREAT! !  Ithere  are  no  errors.') 
error_no=0; 
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errmx=errors; 
rowerrs=sum(  errors'); 
else 

disp('WARNING!:Errors  were  detected!') 

dispC) 

if  n==k 

dispCWARNINGl:  Since  n=k,Conv.  coding  has  no  error  correcting  possible') 

disp(") 

end 

disp(['For  the  given  input  parameters:n=',int2str(n),'and  k=',int2str(k),',the  convolu¬ 
tional  code  is  capable']) 

disp(['of  correcting  ',int2str(symcorr), 'errors.']) 
dispC) 

%Check  to  see  if  xbit  and  ybit  are  the  same.  If  not,  display  error  message 
if  xbit==ybit; 

dispCGREATl! Ithere  are  no  bit  errors.') 

bit_error_no=0; 

bit_errmx=bit_errors; 

b  it_rowerrs=sum(bit_errors') ; 

else 

%disp('WARNING!:Errors  were  detected!') 

dispC) 

if  n==k 

disp('WARNING!:  Since  n=k,Conv.  coding  has  no  error  correcting  possible') 
dispC) 
end 
end 

%CONV  code  was  able  to  correct  all  errors 
if  errtot==0 

Pre_RS_error_matrix=errors; 

disp('EXCEEEENT;  Convolutional  code  corrected  all  detected  errors!') 
disp( ['Originally  the  error  total  was:',int2str(Error_number)]) 
dispC) 
error_no=0; 

[rx  cx]=size(x); 

errmx=zeros(rx,cx); 

rowerrs=sum(errmx'); 

%CONV  code  was  able  to  correct  some  errors  but  not  all  of  them 
elseif  errtoKErrornumber 
Pre_RS_error_matrix=errors; 

Post_RS_error_matrix=RSerrs; 

errmx=RSerrs; 

rowerrs=sum(errmx'); 

dispCOOOPS:  Convolutional  code  corrected  some  detected  errors,  but  not  all.') 

disp( ['Originally  the  error  total  was  :  ',int2str(Error_number)]) 

dispC) 
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disp([' After  Convolutioanl  decoding  ,  the  error  number  was  reduced 
to:',int2str(RSerrtot)]) 
disp(") 

error_no=RSerrtot; 

disp(['The  total  number  of  correct  symbols  are:',mt2str((size(y,l)*size(y,2))-RSerrtot)]) 
disp(") 

%disp('The  error  number  distribution  per  block  number  is  :') 

%disp(errperblk) 

%figure(pic+3) 

%bar((l  :blknos),errblksum) 

%axis([0.5(blknos+0.5)0(max(errblksum)+l)]) 

%title(['Simulation#',int2str(pic),':Error  Distribution  Per  Message  Block  (Error 
count=',int2str(error_no),')']) 

%xlabel( ['Message  Block  Number(block  size:', int2str(blklgth), ’symbols)']) 

%conv  code  did  not  correct  any  errors 
elseif  errtot==Error_number 
Error_matrix=errors; 
errmx=errors; 
rowerrs=sum(  errors'); 

disp('OOOPS!:The  convolutional  code  did  not  correct  any  errors.') 

disp('Perhaps  a  more  powerful  code  is  required.') 

dispC) 

disp(['The  total  number  of  error  occurrences  is:',int2str(Error_number)]) 
disp(") 

error_no=errtot; 

%disp('The  error  number  distribution  per  block  number  is  :') 

%disp(errperblk) 

%figure(pic+4) 

%bar((l  :blknos),errblksum) 

%axis([0.5  (blknos+.5)  0  (max(errblksum)+l)]) 

%title(['Simulation#',int2str(pic),':Error  Distribution  Per  Message  Block.  (Error 
count=',int2str(error_no),')']) 

%xlabel( ['Message  Block  Number  (block  size:', int2str(blklgth)', symbols)']) 

% 

end 

end 

end 

end 

end 

dispC _ '); 

elseif  ecc==2; 
k=l; 

for  i=l:r_rm; 

k=k+prod(l  :m_rm)/prod(l  :m_rm-i)/prod(l  :i); 
end 
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disp(['The  rate  of  this  code  is  ',int2str(k),'  /  ',int2str(2^m_rm),'  = 
',num2str(k/(2^m_rm))]) 
y=y(l,[l;size(x,2)]); 
ybit=ybit(l  ,[1  :size(xbit,2)]); 

%Compare  inputs  x  and  y  and  generate  error  matrix,  "errors" 
errors=(x~=y); 

%Compare  inputs  xbit  and  ybit  and  generate  BIT  error  matrix,  "bit_errors" 
bit_errors=(xbit~=ybit); 

%Find  the  error  distribution  vs.  OFDM  frequencies 
freqerrs=sum(  errors); 

Error_number=sum(sum(  errors)); 
bit_error_total=sum(sum(bit_errors)); 
e=2^(m_rm-r_rm- 1 )- 1 ; 

%Reed-Muller  symbol  correction  for  k  symbols 
if  m_rm-r_rm<=l 

disp('WARNING!  This  code  does  NOT  have  error  correction  capability! !!') 
elseif  m_rm-r_rm>l 

disp(['This  code  can  correct  up  to  ',int2str(e),'  error(s)  in  a  blocklength  of 
',int2str(2^m_rm),'  bits']) 
end 

errtrans=errors'; 

%Reshape  the  error  matrix  as  a  vector  of  errors 
errvect=errtrans(:)'; 
blkrem=rem(length(errvect),k); 
if  blkrem~=0; 
zeropad=zeros(k-blkrem); 
errvectpad=[errvect  zeropad(l,:)]; 
elseif  blkrem==0; 
errvectpad=errvect; 
end 

blknos=length(errvectpad)/k; 
errcorct=[]; 
errblksum=[]; 
for  lp=l:blknos; 

errblk=errvectpad(((k*(lp-l))+l):(k*lp)); 

errblklgth=length(errblk); 

if  sum(errblk)<=e; 

noerr=zeros(errblklgth); 

errblk=noerr(l,:); 

elseif  sum(errblk)>e; 

errblk=errblk; 

end 

errcorct=[errcorct  errblk]; 
errblksum=[errblksum  sum(errblk)]; 
end 

newerrvect=errcorct(  1  :length(errvect)); 
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errtot=sum(newerrvect) ; 

RMerrs=(reshape(newerrvect,size(errors,2),size(errors,l)))'; 

%Find  the  error  distribution  vs.  OFDM  Frequencies 
freqerrs=sum(RMerrs); 
errindex=(find(RMerrs))'; 

RMerrtot=sum(errblksum); 

RMerrdif=Error_number-RMerrtot; 
errperblk=[(l  :blknos);errblksum]; 

%Check  to  see  if  x  and  y  are  the  same.  If  not,  display  error  message 
if  x==y; 

disp('GREAT! !  Ithere  are  no  errors.') 

error_no=0; 

errmx=errors; 

rowerrs=sum(  errors'); 

else 

disp('WARNING!:Errors  were  detected!') 
dispC) 

%Check  to  see  if  xbit  and  ybit  are  the  same.  If  not,  display  error  message 
if  xbit==ybit; 

dispCGREATl! Ithere  are  no  bit  errors.') 

bit_error_no=0; 

bit_errmx=bit_errors; 

b  itro  werrs=sum(bit_errors') ; 

else 

%disp('WARNING!:Errors  were  detected!') 
dispC) 

%RM  code  was  able  to  correct  all  errors 
%if  Error_number==0 
if  errtot==0 

disp('EXCEEEENT:  The  Reed-Muller  code  corrected  all  detected  errors!') 

disp( ['Originally  the  error  total  was:',int2str(Error_number)]) 

dispC) 

error_no=0; 

[rx  cx]=size(x); 

errmx=zeros(rx,cx); 

rowerrs=sum(errmx'); 

%RM  code  was  able  to  correct  some  errors  but  not  all  of  them 
elseif  errtoKErrornumber 
Pre_RM_error_matrix=errors ; 

Post_RM_error_matrix=RMerrs; 

errmx=RMerrs; 

rowerrs=sum(errmx'); 

disp('The  Reed-Muller  code  corrected  some  detected  errors,  but  not  all.') 
disp( ['Originally  the  error  total  was  :  ',int2str(Error_number)]) 
dispC) 

disp([' After  R-M  decoding  ,  the  error  number  was  reduced  to:',int2str(RMerrtot)]) 
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disp(") 

error_no=RMerrtot; 

disp(['The  total  number  of  correct  symbols  are:',mt2str((size(y,l)*size(y,2))- 
RMerrtot)]) 
disp(") 

%disp('The  error  number  distribution  per  block  number  is  :') 
%disp(errperblk);  %degisme  degisme  ;  koyulacak 
%figure(pic+3) 

%bar((l  :blknos),errblksum) 

%axis([0.5(blknos+0.5)0(max(errblksum)+l)]) 

%title(['Simulation#',mt2str(pic),':Error  Distribution  Per  Message  Block  (Error 
count=',mt2str(error_no),')']) 

%xlabel( ['Message  Block  Number(block  size:', mt2str(blklgth), ’symbols)']) 
%RM  code  did  not  correct  any  errors 
elseif  errtot==Error_number 
Error_matrix=errors; 
errmx=errors; 
rowerrs=sum(  errors'); 

disp('The  Reed-Muller  code  did  not  correct  any  errors.') 
disp('Perhaps  a  more  powerful  R-M  code  is  required.') 
dispC) 

disp(['The  total  number  of  error  occurrences  is:',mt2str(Error_number)]) 
disp(") 

error_no=errtot; 

%disp('The  error  number  distribution  per  block  number  is  :') 
%disp(errperblk) 

%figure(pic+4) 

%bar((l  :blknos),errblksum) 

%axis([0.5  (blknos+.5)  0  (max(errblksum)+l)]) 

%title(['Simulation#',mt2str(pic),':Error  Distribution  Per  Message  Block.  (Error 
count=',mt2str(error_no),')']) 

%xlabel( ['Message  Block  Number  (block  size:', mt2str(blklgth)', symbols)']) 

dispC _ '); 

end 

end 

end 

end 

% - 
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CHUHF 


%CHUHF 

%  TITLE  :  UHF  CHANNEL  MODEL  (MULTIPATH  CHANNEL) 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  y=ohuhf(s,x,loss,dly,dop,N,freqspaoe) 

o=10A(-loss./20); 

deltat=l/(N*freqspaoe); 

d=(dly.  *  .00000 1  )./deltat; 

e^dop./freqspaee; 

[L,Nt]=size(x); 

D=length(d); 

x=x.'; 

x=x(:).'; 

%D  path  with  delays  from  d.  (Uses  maero  dline.m) 
xd=dline(x,d); 

[rr,ee]=size(xd); 

x=xd(l,:); 

%  Offsets  direet  path  by  .7  of  max  doppler  fireq.  (uses  maero  ofst.m) 
xo=ofst(.7*e(l),N,x); 

%  First  path  with  no  fading,  (uses  maero  ray  dop.m) 
for  1=1  :D 

a=ray_dop(s,oo,N,e(l)); 

xd(l,:)=a.*xd(l,:); 

end 

%Sums  the  fading  paths 
y=e*xd; 

%Adds  in  the  First  path  without  fading  for  the  GSM-Rieean. 

% 

y=y+xo; 

y=y(l:L*Nt); 

y=reshape(y,Nt,L).'; 

% - 
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CMV2FA 


%CMV2FA 

%  TITLE  :  COMPLEX  EREQUENCY  ARRAY  GENERATOR 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Dave  Roderiek,  Naval  Postgraduate  Sehool 

%  Modified  By  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  X=omv2fa(N,M) 

[m  n]=size(M); 
if  rem(n,2)==0; 

M=M; 

else 

M=[zeros(m,l)  M]; 
end 

[m  n]=size(M); 

K=round(n/2); 

%Generate  a  matrix  of  zeros  with  m  row  and  N  eolumns. 
X=zeros(m,N); 

X(:,1:K)=M(:,K+1:2*K); 

X(:,N-K+1:N)=M(:,1:K); 

% - 
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CNV  ENCD 


%CNV_ENCD 

%  TITLE  :  CONVOLUTIONAL  ENCODING 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  oe_output=onv_enod(oe_g,oe_kO,oe_input) 

%  onv_enod(oe_g,oe_kO,oe_input) 

%  determines  the  output  sequenee  of  a  binary  eonvolutional  eneoder 
%  ee  g  is  the  generator  matrix  of  the  convolutional  code 
%  with  ce  nO  rows  and  ce_l*ce_kO  columns.  Its  rows  are  ce_gl,ce_g2,....ce_gn. 
%  ce  kO  is  the  number  of  bits  entering  the  encoder  at  each  clock  cycle. 

%  check  to  see  if  extra  zero  padding  is  necessary 
if  rem(length(ce_input),ce_kO)>0 

ce_input=[ce_input,zeros(size(l:ce_kO-rem(length(ce_input),ce_kO)))]; 

end 

ce_n=length(ce_input)/ce_kO; 

%check  the  size  of  matrix  ee  g 
if  rem(size(ce_g,2),ce_k0)>0 
error('Error,  ee  g  is  not  of  the  right  size.') 
end 

%  determine  ce_l  and  ce  nO 
ce_l=size(ce_g,2)/ce_k0; 

%disp(['The  value  of  ce_l  is:',mt2str(ce_l)]); 
ce_nO=size(ce_g,  1 ); 

%disp(' ') 

%disp(['The  value  of  ce  nO  is:',mt2str(ce_n0)]); 

%add  extra  zeros 

ce_u=[zeros(size(l;(ce_l-l)*ce_kO)),ce_input,zeros(size(l:(ce_l-l)*ce_kO))]; 

%generate  ce_uu,  a  matrix  whose  column  are  the  contents  of 

%conv.  encoder  at  various  cycles. 

ce_ul=ce_u(ce_l*ce_kO;-l ;  1); 

for  ce_i=l:ce_n+ce_l-2 

ce_u  1 =[ce_u  1  ,ce_u((ce_i+ce_l)  *  cekO :  - 1 :  ce_i*  ce_kO+ 1 )] ; 
end 

ce_uu=reshape(ce_ul,ce_l*ce_kO,ce_n+ce_l-l); 

%determine  the  ce_output 

ce_output=reshape(rem(ce_g*ce_uu,2),l,ce_n0*(ce_l+ce_n-l)); 

% - 


84 


COFDMSIM 


%COFDMSIM 

%  TITLE  ;  SIMULATION  OF  COFDM 

% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Dave  Roderiek,  Naval  Postgraduate  Sehool 

%  Modified  By  :  TEZEREN,Serdar  Umit,FTJG, TURKISH  NAVY 

% - 

ele,elose  all;elear  all; 

dispC _ '); 

disp('This  bateh  m-file  runs  COFDM  simulations  using  different  ehannel  models.') 
fort=input('To  run  the  frequeney  version,  enter  l(one),  To  run  the  time  version,  enter 
O(zero),  or  to  run  both  enter  2(two):'); 

freqno=input('Enter  the  #  of  OFDM  frequeneies  (note  :  must  be  even):'); 
%N=input('Enter  the  number  of  FFT  points  (Note  :  This  number  must  be  larger  than  #  of 
OFDM  frequeneies):'); 

N=64; 

ehnmdl=input('Do  you  want  to  run  ehannel  model  0,  1,  2,  3  or  4?  (Enter  0,1, 2,3  or  4):'); 
if  ohnmdl==0 

dispCCode  eheek  simulation  performed.'); 

sigs=0; 

loss=0; 

dop=0; 

dly=0; 

elseif  chnmdl==1 

disp('AWGN  Channel  simulation  performed.'); 

sigs=input('Enter  the  sigma  noise  parameter  range  or  single  value.  (Ex  lin- 

spaee(0,0.02,20)or  .003):'); 

loss=0; 

dop=0; 

dly=0; 

elseif  ehnmdl==2 

disp('Multipath  Channel  simulation  performed.'); 
sigs=0; 

pthno=input('Do  you  want  to  run  linkl  ,link2,  link3  or  a  eustom  link  ?  (Enter  1,2,3  or  4 
for  eustom):'); 

% 

%Eink  parameters 

% 

if  pthno==3 
%my  link  3 
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loss=[0, 2.17, 4.34, 6.51, 8.69, 10.86, 13.03, 15.20, 17.37, 19.54, 21.71, 23. 89, 26.06, 28.23, 30.4, 

32.57,34.74,36.92]; 

dop=[15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15]; 

dly=[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0. 

85]; 

elseif  pthno==2 
%my  link  2 

loss=[0,2.17,4.34,6.51,8.69,10.86,13.03, 15.20,17.37,19.54,21.71,23. 89,26.06,28.23,30.4, 

32.57,34.74,36.92]; 

dop=[10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10]; 

dly=[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0. 

85]; 

elseif  pthno==l 
%my  link  1 

loss=[0,2.17,4.34,6.51,8.69,10.86,13.03, 15.20,17.37,19.54,21.71,23. 89,26.06,28.23,30.4, 

32.57,34.74,36.92]; 

dop  1^5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5], 

dly=[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0. 

85]; 

elseif  pthno==4 

disp('Custom  link  simulation...') 
loss=input('Enter  the  path  loss  in  db  (Ex 

[0,2.17,4.34,6.51,8.69,10.86,13.03,15.20,17.37,19.54,21.71,23.89,26.06,28.23,30.4,32.57 

,34.74,36.92]):'); 

dop=input('Enter  the  doppler  frequeney  in  Hertz  (Ex 
[5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5 ]).  ), 
dly=input('Enter  the  time  delays  of  the  multipaths  in  mierosees  (Ex 
[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0.85]):' 
); 

end 

elseif  ehnmdl==3 

disp('AWGN+Multipath  Channel  simulation  performed'); 
sigs=input('Enter  the  sigma  noise  parameter  range  or  single  value.  (Ex  lin- 
spaee(0, 0.02,20)  or  .003):'); 

pthno=input('Do  you  want  to  run  linkl,  lmk2,  lmk3  or  a  eustom  link  ?(Enter  1,2,3  or  4 
for  eustom):'); 

% 

%Eink  parameters 
% 

if  pthno==3 
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%my  link  3 


loss=[0, 2.17, 4.34, 6.51, 8.69, 10.86, 13.03, 15.20, 17.37, 19.54, 21.71, 23. 89, 26.06, 28.23, 30.4, 

32.57,34.74,36.92]; 

dop=[15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15]; 


dly=[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0. 

85]; 

elseif  pthno==2 
%my  link  2 


loss=[0,2.17,4.34,6.51,8.69,10.86,13.03, 15.20,17.37,19.54,21.71,23. 89,26.06,28.23,30.4, 

32.57,34.74,36.92]; 

dop=[10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10]; 


dly=[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0. 

85]; 

elseif  pthno==l 
%my  link  1 


loss=[0,2.17,4.34,6.51,8.69,10.86,13.03, 15.20,17.37,19.54,21.71,23. 89,26.06,28.23,30.4, 

32.57,34.74,36.92]; 

dop  1^5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5], 


dly=[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0. 

85]; 

elseif  pthno==4 

disp('Custom  link  simulation...') 
loss=input('Enter  the  path  loss  in  db  (Ex 

[0,2.17,4.34,6.51,8.69,10.86,13.03,15.20,17.37,19.54,21.71,23.89,26.06,28.23,30.4,32.57 

,34.74,36.92]):'); 

dop=input('Enter  the  doppler  frequency  in  Hertz  (Ex 
[5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5 ]).  ), 
dly=input('Enter  the  time  delays  of  the  multipaths  in  microsecs  (Ex 
[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0.85]):' 


); 

end 

elseif  chnmdl==4; 

sigs=input('Enter  the  sigma  noise  parameter  range  or  single  value.  (Ex  lin- 
space(0,0.02,20)or  .003):'); 

disp('Mobile  channel  simulation') 

%sigs=0; 

loss=0; 

dop=0; 

dly=0; 
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end 

%%%%allcase=input('Simulate  all  interleaver  eases  (yes)  or  speeifie  ones(no)? 

(l=yes,0=no):'); 

allease=0; 

if  allease==l 

disp('All  eases, (0-8), will  be  tested.'); 

eases=[0:8]; 

elseif  allease==0 

%%%%eases=input('Enter  speeifie  ease  numbers  from  (0  to  8)(Ex  [0  4  5  8]):'); 

oases=0; 

end 

if  fort~=2 

if  length(eases)'-'=l 

easey_n=input('Do  you  want  to  find  optimal  interleaver  ease(s)  ?  (l=yes,  0=no):'); 
else 

easey_n=0; 

end 

end 

totsym=input('Enter  the  total  minimum  number  of  symbols  to  simulate  (Ex  10000):'); 

rowno=eeil(totsym/freqno); 

if  totsym~=(rowno*freqno) 

disp(['Note:Based  on  the  parameters  thus  far,  the  aetual  total  number  of  symbol  to  be 

simulated  will  be  :',int2str(rowno*freqno)]); 

end 

%pry_n=input('Eor  the  interleaver,  do  you  want  to  ealeulate  all  possible  intermediate  ma¬ 
trix  dimension  pairs?(l=yes,0=no):'); 
pry_n=l; 
pairl=l; 

pair2=ro  wno  *  fireqno ; 
if  pry_n==l 

Intrlvr_pairs=intlvprs(ro  wno ,  freqno) ; 
intlvrprs=Intrlvr_pairs; 

%ee=round(size(intlvrprs,l)/3); 

%int( : ,  1  )=intlvrprs(  [  1 :  ee] ,  1 ) ; 

%mt( :  ,2)=mtlvrprs(  [  1 :  ee]  ,2) ; 

%int( :  ,3  )=intlvrprs(  [ee+ 1 : 2  *  ee] ,  1 ) ; 

%mt( :  ,4)=mtlvrprs(  [ee+ 1 :2  *  ee]  ,2); 

%int( :  ,5)=intlvrprs(  [2  *  ee+ 1 :  end] ,  1 ); 

%int( :  ,6)=intlvrprs(  [2  *  ee+ 1 :  end]  ,2); 
dispC) 

disp('Eor  these  input  parameters,  all  possible  aeeeptable  interleaver  dimension  pairs  are: ') 
disp(Intrlvr_pairs) 

%disp(int) 

end 

pairs=mput(['Desired  interleaver  pair?  (Ex  [row  #  eol  #]  =  [20  50]  ):']); 
rmtlv=pairs(l); 
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cintlv=pairs(2); 

%mary=input('Enter  the  number  of  M-ary  bits,  q  (i.e.  for  256-ary,  q=8):'); 
mary=l; 

%nary=input('Enter  the  number  of  N-ary  bits,q(i.e.  for  16-ary,  q=4):'); 
nary =2; 

freqspaoe=round(  1 6600000/freqno); 

%Ng=input('Enter  the  guard  interval  length  (Number  of  sample  points):'); 

Ng=16; 

ecc=input('Enter  the  type  of  error  correetion  eoding  ?  (0=No  EEC  code,  l=Convolutional 
Code,  2=Reed-Muller  Code):'); 
if  ecc==l 

code=input('Enter  n,  k  and  error  correction  block  length  (Ex:  [240  200  240]):'); 

n=code(l); 

k=code(2); 

blklgth=code(3); 

r_rm=0;m_rm=0; 

elseif  ecc==2; 

pair_rm=input('Enter  the  r,m  pair  for  the  Reed-Muller  code  (Ex:  [2  4]):'); 
r_rm=pair_rm(  1 , 1  );m_rm=pair_rm(  1 ,2); 
if  r_rm>m_rm 

display('ERROR! ! !  m  must  be  equal  to  or  larger  than  r,  and  both  must  be  integer  larger 
than  zero!') 

pair_rm=input('Enter  the  r,m  pair  for  the  Reed-Muller  code  (Ex:  [2  4]):'); 

r_rm=pair_rm(  1 , 1  );m_rm=pair_rm(  1 ,2); 

end 

n=0;k=0;blklgth=0; 

elseif  ecc==0 

n=freqno; 

k=freqno; 

blklgth=freqno; 

r_rm=0;m_rm=0; 

end 

%svals=input('Enter  specific  seed  values,  or  0  for  a  random  seed  (ex  [103  22,  60]  or 
[0]);’); 
svals=33; 

%picy_n=input('Do  you  want  signal  plots?  (l=yes,  0=no):'); 

picy_n=l; 

if  picy_n==l 

%wait=input('How  many  seconds  of  delay  between  pictures?'); 
wait=3; 

wait=round(wait) ; 
elseif  picy_n==0 
wait=0; 
end 

%prnty_n=input('Do  you  want  print  outs?  (l=yes,  0=no):'); 
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pmty_n=0; 

pic=0; 

svect=[]; 

for  run=l  :length(svals); 
errvect=[]; 
errcase=[]; 
errtot=[]; 
if  min(svals)==0 

rand('seed',sum(100*clock)); 
s=round(abs(rand(  l)*pi*10*(pic+ 1  )*run)); 
elseif  min(svals)~=0 
s=svals(run); 
end 

svect=[svect,s]; 
for  l=l;length(cases); 

dispC _ ') 

disp(['Run#:',int2str(run)]); 
if  chnmdl==0 

disp( ['Channel  Model  0  :  Code  Cheek  Simulation  Performed']) 
elseif  chnmdl==l 

disp(  ['Channel  Model  1  :  AWGN  Channel  with  Performed']) 
elseif  ehnmdl==2 

disp( ['Channel  Model  2  Link  ',int2str(pthno),' :  Multipath  Channel  Performed']) 
elseif  ehnmdl==3 

disp( ['Channel  Model  3  Link  ',int2str(pthno),' :  AWGN  +  Multipath  Channel  Per¬ 
formed']) 

elseif  ehnmdl==4 

disp( ['Channel  Model  4  :  Severe  Mobile  Outdoor  Channel  Performed']); 
end 

%disp(['Seed=',int2str(s)]); 
disp(  ['Interleaver  ease=',int2str(eases(l))]); 
disp(  ['Interleaver  Pair=[',int2str(pairs), ']']); 
disp(['The  #  of  OFDM  suboarriers=',int2str(freqno)]); 
if  eoo==0; 

disp(['FEC  eode  is  not  used  for  this  run']); 
elseif  eee==I 

disp(['n=',int2str(n),'  k=',int2str(k),'  eonvolutional  eode  is  used  for  this  run']); 
elseif  eee==2 

disp(['R(',int2str(r_rm),',',mt2str(m_rm),')  Reed-Muller  eode  is  used  for  this  run']); 
end 

if  fort<=I 

[errmax, errors, freqerrs,papr]=ehanedl(ehnmdl, wait, prnty_n,piey_n, pie, eases(l),s,freqno,ri 
ntlv,eintlv,N,mary,nary,n,k,blklgth,Ng,sigs,loss,dly,dop,freqspaee,fort,eee,r_rm,m_rm); 
elseif  fort==2 

disp('Frequeney  differential  eneoding/deeoding  simulation...') 

dispC _ ') 
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[errmax, errors, freqerrs,papr]=chancdl(chnmdl, wait, prnty_n,picy_n, pic, cases(l),s,freqno,ri 
ntlv,cintlv,N,mary,nary,n,k,blklgth,Ng,sigs,loss,dly,dop,freqspace,l,ecc,r_rm,m_rm); 

disp('Time  differential  encoding/decoding  simulation....') 

dispC _ ') 

[errmax,errors,freqerrs,papr]=cliancdl(chnmdl,wait,prnty_n,picy_n,pic+12,cases(l),s,freq 

no,rintlv,cintlv,N,mary,nary,n,k,blklgth,Ng,sigs,loss,dly,dop,freqspace,0,ecc,r_rm,m_rm) 

end 

errtot=[errtot  sum( errors)]; 
errvect=[errvect  errtot]; 
errcase=[errcase  sum( errmax)]; 
end 

if  fort~=2 

casearry=[cases;errcase]; 
if  casey_n==l 
figure(pic+14) 
bar(cases,errcase) 
grid 

orient  tall 
if  fort==l 

title([int2str(pic),':Maximum  Error  Total  Vs.  Interleaver  CASE  Number 
(Ereq.Diff.  Enc.)  (OEDM  Ereq.#=',int2str(freqno),')']) 
elseif  fort==0 

title([int2str(pic),':Maximum  Error  Total  Vs.  Interleaver  CASE  Number  (Time 
Diff.  Enc.)  (OEDM  Ereq.#=',int2str(freqno),')']) 
end 

xlabel(['CDE  Interleaver  CASE  Number']) 

ylabel(['Maximum  Error  Count  Eor  Any  Symbol  Row  (Seed=',int2str(s),')']) 
axis([-.5  8.5  0  (max(errcase)+l)]) 
if  pmty_n==l; 
print 

pause(lO) 

end 

pause(wait); 

figure(pic+15) 

bar(cases,errtot) 

grid 

orient  tall 

title([int2str(pic),':Error  Totals  Vs.  Interleaver  CASE  Number']) 
xlabel(['CDL  Interleaver  CASE  Number']) 

ylabel(['Sigma:(',num2str(min(sigs)),'-',num2str(max(sigs)),')  Error  Total']) 
axis([-.5  8.5  (min(errtot)-l)  (max(errtot)+l)]) 
if  pmty_n==l; 
print 
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pause(lO) 

end 

pause(wait); 

end 


pie=pie+l; 

end 

end 

disp(") 

disp('Channel  model  bateh  run  is  finished!') 
%Seed=sveet 


%■ 
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CVDD 


%CVDD 

%  TITLE  :  CONTINUOUS  VARIABLE  DIGITAL  DELAY  ELEMENT 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  [y]=cvdd(x, alpha) 
if  ((nargin~=2)|(nargout~=l)) 
error('ERROR:usage:y=y=ovdd(x, alpha);'); 
return; 
end 

if  (size(x)~=size(alpha)) 

error('ERROR:x  and  alpha  must  be  the  same  size'); 

return; 

end 

if  (abs(alpha)>0.5) 

error('ERROR: alpha  must  be  within  -0.5  and  0.5'); 

return; 

end 

% - 

%  Initialization 

% - 

%  Initialize  EIR  filter  eoeffieients  are  in  [1]  (0,0.328  pass  band) 

C0=[-0.013824  0.054062  -0.157959  0.616394  0.616394  -0.157959  0.054062  -0.013824]; 
Cl=[0.003143  -0.019287  0.1008  -1.226364  1.226364  -0.1008  0.019287  -0.003143]; 
C2=[0.055298  -0.216248  0.631836  -0.465576  -0.465576  0.631836  -0.216248  0.055298]; 
C3=[-0.012573  0.077148  -0.403198  0.905457  -0.905457  0.403198  -0.077148  0.012573]; 

% - 

%  4  parallel  EIR  and  add  together  based  on  [1] 

% - 

y0=fdter(C0,[l],x); 

yl=fdter(Cl,[l],x); 

y2=fdter(C2,[l],x); 

y3=fdter(C3,[l],x); 

% 

y=alpha.*y3; 

y=alpha.*(y+y2); 

y=alpha.*(y+yl); 

y=y+y0; 

% - 
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DECDRCDL 


%DECDRCDL 
%  Title  :  COFDM  DECODER 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Dave  Roderiek,  Naval  Postgraduate  Sehool 

%  Modified  By  :  TEZEREN,Serdar  Umit,ETJG, TURKISH  NAVY 

% - 

fune- 

tion[outmsg,deeoder_output_bit,random_msg,random_bit,M, MM, binary  _value]=deedred 
l(piey_n, pie, dease,K, Fa, nsymno,freqno,rdintlv,edintlv,mary, nary, fort, B_random,eee,r_rm 
,m_rm) 

%to  generate  BER  vs  Eb/No. 

%to  provide  bit  outputs 
M=fa2oma(K,Fa); 

Cmp  lx_mod_vals=M ; 

% 

naryp=nary; 

[s,MM]=dfdodrft(naryp,nary,M,fort); 

[E,co]=size(s); 
strans=s'; 
sveot=strans(:).'; 
oorrs=sveot(l  insymno); 

% 

nsymno; 

Br=bm(mary,mb(nary,oorrs)); 
lengthBr=length(Br); 
rmndr=rem(length(Br),freqno); 
if  rmndr==0; 

Br=Br; 

elseif  rmndr~=0; 

Br=Br(  1  ;(lengthBr-rmndr)); 
end 

re  vd=(reshape(Br,  freqno ,  length(Br)/ fireqno))' ; 

Revd_Intlv_Ary=revd; 

% 

[Br  Be]=size(revd); 

SYNC=[]; 

sr=revd'; 

si=sr(:)'; 

sd=edldlv(rdintlv,edintlv,dease,si,SYNC); 
binary_value=mb(mary,sd); 
if  eee~=2; 
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viter_G=[l  0  1  1  0  1  1;1  1  1  1  0  0  1]; 
viter_k=l; 

[viterbi_output,survivor_sta,cumul_metrix]=viterbi(viter_G,viter_k,binary_value); 

decoder_output=viterbi_output; 

%mary_dec=bm(mary,decoder_output); 

%decoder_output_bit=decoder_output; 

%outmsg=reshape(sd,Bc,Br)'; 

% 

%random_bit=B_random; 

%random_msg=bm(mary,random_bit); 

%[Brow  Bcol]=size(random_msg); 

% 

%outmsg=reshape(mary_dec,Bcol,Brow)';bu  zaten  yuzdeliydi  kaldirma 

%Sink_Msg=outmsg; 

elseif  ecc==2; 

[rmdecoder_output]=rmdecoder(binary_value,r_rm,m_rm); 

decoder_output=rmdecoder_output; 

end 

mary_dec=bm(mary,decoder_output); 

decoder_output_bit=decoder_output; 

%outmsg=reshape(sd,Bc,Br)';  bu  zaten  yuzdeliydi  kaldirma 
% 

random_bit=B_random; 

random_msg=bm(mary,random_bit); 

[Brow  Bcol]=size(random_msg); 

% 

outmsg=reshape(mary_dee,size(mary_dec,l)*size(mary_dee,2),l)'; 

Sink_Msg=outmsg; 

% - 
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DECI2BIN 


%DECI2BIN 

%  TITLE  :  DECIMAL  TO  BINARY  CONVERSION 

% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 

%  Author  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  y=deoi2bin(x,l) 

y=zeros(l,l); 

vi=l; 

while  x>=0  &  vi<=l 

y(vi)=rem(x,2); 

x=(x-y(vi))/2; 

vi=vi+l; 

end 

y=y(l:-i:i); 

% - 
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DFDCDRFT 


%DFDCDRFT 

%  TITLE  :  COMPLEX  NUMBER  DEMODULATOR  &  EREQUENCY  /  TIME 

DIEEERENTIAL  DECODING 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  [s,M]=dfdodrft(qp,q,MD,fort) 
if  fort==0  %Time  Differential  deeoding 
MD=MD’; 

[m  n]=size(MD); 

%  Perform  a  looping  routine  to  find  the  phase  differenees  between  adjaeent  values  in  the 
%  array, MD, and  put  these  ealeulated  values  into  array ,M. 
for  1=1  :m 
for  j=l;n-l 

M(l,j)=MD(l,j+l)*eonj-(MD(l,j)); 

end 

end 

%Transpose  the  array  baek  to  its  original  form 
M=M’; 

%  Caleulate  the  number  of  M-ary  symbols  based  upon  the  exponent  qp,then  use  this 
number 

%  to  find  the  number  of  equally  spaeed  phases  in  a  unit  eirele. 

N=2^qp; 

dph=2*pi/N; 

%  Divide  the  phase  arguments  of  elements  in  M,  by  the  equal  phases  generated  by  dph. 
phn=angle(M)  ./dph; 

%  Calculate  the  phase  sector  number  by  finding  the  remainders. 
s=rem(round(phn)+N,N); 

elseif  fort==l  %  Erequency  Differential  decoding 
%  Transpose  the  modulation  array,  and  find  the  dimensions 
[m,n]=size(MD); 

MD=MD(:,2:n); 

[m  n]=size(MD); 

%  Perform  a  looping  routine  to  find  the  phase  differences  between 
%  adjacent  values  in  the  array,  MD,  and  put  these  calculated  values  into  array ,M. 
for  1=1  :m 
for  j=l;n-l 

M(l,j)=MD(l,j+l)*conj(MD(l,j)); 

end 

end 
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N=2^qp; 

dph=2*pi/N; 

%  Calculate  the  phase  seetor  number  by  finding  the  remainders. 

phn=angle(M)  ./dph; 

s=rem(round(phn)+N,N); 

end 

% - 
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DIFCDRFT 


%DIFCDRFT 

%  TITLE  :  COMPLEX  NUMBER  MODULATOR  &  EREQUENCY  /  TIME 
DIEEERENTIAL  ENCODING 


% - 

%  Thesis  Advisor 
%  Author 
%  Modified  by 
% - 


:  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 
;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 


funetion  MD=difodrft(q,m,fort) 
if  fort==0  %Time  differential  eneoding 
%  M-ary  alphabet  size 
N=2^q; 

%  Determine  the  number  of  equal  phases  based  upon  the  m-ary  symbol  length 
dph=2*pi/N; 

%  find  the  size  of  the  input  symbol  matrix  (#  of  row  &  #  of  eolumns) 

[rr  n]=size(m); 

%  Perform  the  time  differential  eneoding  of  phase  values  by  eumulative  summing  matrix, 
%  m,  down  one  eolumn  at  a  time  aeross  the  entire  matrix.  This  funetion  generates  a  ma¬ 
trix. 

for  k=l;n 

md=oumsum(m( :  ,k)) ; 

%  Generate  the  eomplex  numbers  with  eorrespondiing  phase  values. 

MD(:,k)=exp(i*dph.*md); 

end 

%  Injeet  the  referenee  row  of  ones  (zero  phase)  at  top  of  output  matrix  for 
%  differential  eneoding  synehronization 
MD=[ones(l,n);  MD]; 

elseif  fort==l  %  Erequeney  Differential  eneoding 

%  M-ary  alphabet  size 

N=2^q;dph=2*pi/N; 

%  Eind  the  size  of  the  input  symbol  matrix  (#  of  row  &  #  of  eolumns) 

[rr  n]=size(m); 
md=oumsum(m') ; 
md=md'; 

%  Generate  the  eomplex  numbers  with  eorresponding  phase  values. 

MD=exp(i*dph.*md); 

%  Injeet  the  referenee  row  of  ones  (zero  phase)  at  top  of  output  matrix  for 
%  differential  eneoding  synehronization. 

MD=[ones(rr,2)  MD]; 
end 

% - 
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DLINE 


%DLINE 

%  TITLE  ;  UHF  CHANNEL  DELAY  LINE  GENERATOR 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  xd=dline(x,d) 
x=x.'; 

dmax=max(d); 
dmin=min(d); 
nmin=floor(dmin); 
nmax=oeil(dmax); 
x=[x;zeros(nmax+3 , 1 )] ; 

N=length(x); 

Nd=length(d); 

% 

for  n=l:Nd; 
di=d(n); 

D=floor(di); 

deld=di-D; 

xd(;,n)=ovdd(x,deld-.5); 

xd(:,n)=[zeros(D,l);xd(l;N-D,n)]; 

end 

xd=xd.'; 

[rr,ee]=size(xd); 
xd=xd( :  ,4+nmin:  ee); 

% - 


100 


DOT  PRO 


%DOT_PRO 

%  TITLE  ;  DOT  PRODUCT  OPERATOR 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  TEZEREN,Serdar  Umit,  ETJG,  TURKISH  NAVY 

% - 

funetion  z=dot_pro(X,S) 

a=size(X); 

z=0; 

for  i=l:a(l,2); 

z=xor(z,(X(i)*S(i))); 

end 

% - 
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FA2CMA 


%FA2CMA 

%  TITLE  :  FREQUENCY  ARRAY  TO  COMPEEX  MODUEATION  ARRAY 
CONVERTER 


% - 

%  Thesis  Advisor  ;  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  :  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 
%  Modified  by  :  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  Mm=fa2oma(K,X) 

[m  n]=size(X); 

Mm(:,l:K)=X(:,n-K+l;n); 

Mm(:,K+l:2*K)=X(:,l;K); 

Cmplx_mod_vals=Mm; 

% - 
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GM  PARTI 


%GM_PART1 

%  TITLE  :  RM  ENCODING  MATRIX  EIRST  PART  GENERATOR 


% - 

%  Thesis  Advisor  ;  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  TEZEREN,Serdar  Umit,  ETJG,  TURKISH  NAVY 

% - 

funetion  [bb]=gm_partl(m); 

AA=zeros(2^m,m); 

B=eye(m); 
for  i=2:2^m; 

AA(i,:)=xor(AA(i-l,;),B(m,:)); 
for  j=l:m; 

if  rem(i-l,2^j)==0; 

AA(i,:)=xor(AA(i,;),B(m-j,:)); 

end 

end 

end 

AA=AA'; 
for  i=l:size(AA,l); 
for  i=l:size(AA,2); 
ifAA(i,j)==0; 

AA(i,j)=l; 
elseif  AA(i,j)==l; 

AA(i,j)=0; 

end 

end 

end 

bb=zeros(size(AA,  1 )+ 1  ,size(AA,2)); 
bb(  1 ,  :)=ones(  1  ,size(AA,2)); 
bb([2:end],:)=AA; 

% - 
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INTLVPRS 


%INTLVPRS 

%TITLE  :  INTERMEDIATE  MATRIX  INTERLEAVER  DIMENSION  PAIRS 

% - 

%Thesis  Advisor  ;  Prof  J.  McEachen,  Naval  Postgraduate  School 

%Author  ;  Dave  Roderick,  Naval  Postgraduate  School 

%Modified  By  ;  Tan  Kok  Chye,  Naval  Postgraduate  School 

% - 

function  pairs=intlvprs(n,m) 
prod=n*m; 
multvect=[l]; 
for  i=2;prod; 
remdr=rem(prod,i); 
if  remdr==0 
multvect=[multvect  i]; 
else 

multvect=multvect; 

end 

mult=multvect; 

end 

lngth=length(mult) ; 
nbr=mult(lngth); 
result=[l  nbr]; 
for  i=2:lngth; 

cmtpr=[mult(i)  nbr/mult(i)]; 
result=[result;cmtpr] ; 
end 

pairs=result; 

% - 
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INVERSE 


%TITLE  :  INVERSE  OE  THE  BINARY  ARRAY 


% - 

%  Thesis  Advisor  ;  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  TEZEREN,Serdar  Umit,  ETJG,  TURKISH  NAVY 

% - 

funetion  d=inverse(D) 

e=size(D); 

f=e(I,2); 

d=xor(ones(I,f),D); 

% - 
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ITDA 


%ITDA 

%TITLE  :  IFFT  AND  GI  RFMOVAF  OPERATOR 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  Y=itda(Ng,y) 

[F  Nt]=size(y); 

%  Remove  the  guard  interval  for  ehannel  eompensation,  Ng,  preeursor. 
y=y(:,Ng+l:Nt); 

%  Take  the  FFT  of  array, y 
Y=fft(y.').'; 

% - 
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%JAKES 


%  TITLE  :  RANDOM  VECTOR  GENERATOR  WITH  JAKES  SPECTRUM 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  :  Prof  R.  Cristi,  Naval  Postgraduate  Sehool 

% - 

funetion  g=jakes(Pd,Ps,N); 

%  g=jakes(Pd,Ps,N) 

%  generate  a  random  veetor  of  length  N  with  jakes  speetrum 
%  Pd=doppler  frequeney 
%  Ps=sampling  frequeney  (Ps>10*Pd  required) 

%  N=veetor  length 
%  g=generated  random  veetor 
%  The  veetor  g  has  unit  power,  ie  g'*g*Ps/N=l 
L=round(log2(P  s/(  1 0  *Pd))); 

P0=Ps/(2^L);  %  sampling  frequeney  for  generating  the  random  sequenee 

N0=2*oeil(N/(2^L)); 

Nf=5 12;  %  number  of  frequeney  eomponents. . . 

%  ...  in  speetrum  definition 
Nh=256;  %  length  of  EIR  filter 

%  EIR  Eilter  impulse  response 

w=0.54-0.46*eos(2*pi*(0;Nh-l)/Nh);  %  hamming  window 

fd=Ed/E0;  %  digital  doppler  frequeney 

f=0:(l/Nf):  1-(1/Nf);  %  veetor  of  digital  frequeneies  (l/2=Nyquist  Ereq) 
kd=fioor(fd*Nf);  %  index  for  doppler  frequeney  fd 

H0(1  :kd)=sqrt(l-(f(l  :kd)/fd).^2); 

I=find(H0==0);  H0(I)=0 .0001*  ones(size(I)); 

H=zeros(l,Nf); 

H(l:length(H0))=(l./H0); 

H=H.*exp(-j*2*pi*f*Nh/2); 

hO=real(ifft(H)); 

h=h0(  1 ;  length(w)) .  *  w; 

%  Generate  time  varying  taps  at  low  sampling  freq  Es_oh 
seed=14; 

randn('state',seed); 
a=randn(  1  ,N0+Nh); 
seed=seed+6; 
b=randn(  1  ,N0+Nh); 

%x=randn(  1  ,N0+Nh)+j  *randn(  1  ,N0+Nh); 
x=a+j*b; 

g=filter(h,l,x); 
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g=g(Nh+l  :Nh+NO);  %  steady  state  response 
%  Upsample  to  Fs=(2^L)*F0  in  L  stages 
for  m=l:L 
F0=2*F0; 

g=reshape([g;zeros(size(g))],  1 ,2*length(g)); 
omegad=2*pi*Fd/F0; 

Domega=(pi/2)-omegad; 

omegao=((pi/2)+omegad)/2; 

M=oeil(((8  *pi/Domega)- 1  )/2); 

Nlilt=2*M+l;  %  filter  order  using  hamming  window 

nt=0;Nfilt-l; 

w=0.54-0.46*oos(2*pi*nt/Nfilt);  %  hamming  window 

hm(l  :M)=sin(omegao*(nt(l  :M)-M))./(pi*(nt(l  :M)-M)); 
hm(M+ 1  )=omegao/pi; 

hm(M+2  :N  filt)=sin(omegao  *  (nt(M+2  :Nfilt)-M))  ./(pi  *  (nt(M+2  ;Nfilt)-M)); 

hm=hm(  1  :Nfilt) .  *  w; 

g=filter(hm,l,g); 

g=g(Nfilt:length(g));  %  steady  state 
end 

%  Normalize 
g=g(l:N); 

K=g*g'/N; 

g=g/sqrt(K); 

% - 
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MARYMSG 


%MARYMSG 

%  TITLE  ;  M-ARY  RANDOM  SIGNAL  GENERATOR 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

%  Modified  by  :  TEZEREN,Serdar  Umit,  LTJG,  TURKISH  NAVY 

% - 

funetion  [vmary_ee,random_bit,ene_output]=marymsg(q,n,m,eee,r_rm,m_rm) 
%  if  reedmuller  is  desired 
if  eee==2; 

k=l;  %  #  of  the  rows=k;#  of  the  eolumns=2^m 
for  i=l;r_rm; 

k=k+prod(l  ;m_rm)/prod(l  :m_rm-i)/prod(l  ;i); 
end 

[random_bit]=msg(round(q*n*m*k/(2^m_rm))); 

ene_output=rmeneoder(r_rm,m_rm,random_bit); 

%unless  reedmuller  is  desired 
elseif  eee~=2; 

[random_bit]=msg(n*m*q*.5-6); 
eonv_g=[l  0  1  1  0  1  1;1  1  1  1  00  1]; 
eonv_k0=l; 

ene_output=env_ened(eonv_g,eonv_kO,random_bit); 

end 

deeml=bm(q,ene_output); 
deeiml=deeml(l  ,1  ;(n*m)); 
vmary_ee=(reshape(deeiml,m,n))'; 

% - 
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MB 


%MB 

%  TITLE  :  M-ARY  TO  BINARY  CONVERTER 

% - 

%  Thesis  Advisor  :  Prof  J.  McEachen,  Naval  Postgraduate  Sehool 

%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  [b]=mb(q,m) 

row=size(m,l); 

ool=size(m,2); 

m=reshape(m',  1  ,(row*ool)); 

b0=rem(m,2); 

m=(m-b0)./2; 

B=bO; 
for  j=l;q-l 
bj=rem(m,2); 
m=(m-bj)./2; 

B=[B;bj]; 

end 

b=B(:)'; 

binary=b; 

% - 
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METRIC 


%METRIC 

%  TITLE  :  VITERBI  HARD  DECISION  DECODING  METRIC 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  distanoe=metrio(v_x,v_y) 

if  v_x==v_y 

distanoe=0; 

else 

distanee=l; 

end 

% - 
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MLTPL 


%MLTPL 

%  TITLE  :  COMMON  MULTIPLES 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  [mult]=mltpl(n,m) 
max=n*m; 
multveot=[l]; 
for  i=2:max; 

remdr=rem(max,i); 
if  remdr==0 

multveot=[multveot  i]; 
else 

multveet=multveet; 

end 

mult=multveet; 

end 

% - 
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MSG 


%MSG 

%  TITLE  :  RANDOM  BIT  SEQUENCE  GENERATOR 

% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  TEZEREN,Serdar  Umit,  ETJG,  TURKISH  NAVY 

% - 

funetion  souroe_msg_bit=msg(k) 
seed=33; 

souroe_msg_bit=randint(  1  ,k,[0  1  ],seed); 

%u=zeros(  1  ,k);%u(  1 , 1  )=  1 ; 

% - 
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NXT  STAT 


%NXT_STAT 

%  TITLE  :  NEXT  STATE  IN  VITERBI  AEGORITHM 

% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  [next_state,memory_oontents]=nxt_stat(ourrent_state, input, v_E,v_k) 
binary_state=deoi2bin(ourrent_state,  v_k*  (v_E- 1 )); 
binary_input=deoi2bin(input,v_k); 

next_state_binary=[binary_input,binary_state(l;(v_E-2)*v_k)]; 

next_state=bin2deoi(next_state_binary); 

memory_oontents=[binary_input,binary_state]; 

% - 
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OFST 


%OFST 

%  TITLE  :  CHANNEL  OEESET 

% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  xo=ofst(e,N,x) 

[m  Nt]=size(x); 

xo=x.'; 

x=x.'; 

x=x(:); 

x=x.'; 

Nt=length(x); 

1=1  :Nt; 

ex=x.*exp(i*(2*pi/N)*e.*l); 

xo(:)=x; 

xo=xo.'; 

% - 
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PAPR2 


%PAPR2 

%  TITLE  :  PEAK-TO-AVERAGE  POWER  RATIO,  PEAK  CEIPPING  AND 
HANNING  WINDOWING 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 

%  Author  ;  TEZEREN  Serdar  Umit,  ETJG,  TURKISH  NAVY 

% - 

funetion  [ppp]=papr2(xmtifft,N,freqno,pio); 

window=l;  %0  EORNO  WINDOWING . 1  EOR  WINDOWING...!  EOR  CEIPPING 

if  window==l 
xmtifft=xmtifft  *  5  0 ; 
xmtifft=xmtifft( : ,  [  1 7 : 8  0] ) ; 
sig=abs(xmtifft); 
figure(3  0);plot(sig) 

han=hanning(9)';  %HANNING  WINDOW 
ko=0.1;  %hanning  eoeffieient 
wo=l-ko*han;  %window  construetion 
%size_wo=size(wo) 
ka=0.2; 

wo2=l+ka*han; 

rmspow=sqrt(64); 

cr=1.4; 

A=or*rmspow;  %olipping  level 
sigl=sig'; 
mxm=max(sig  1 ); 
for  i=l:size(sig,l) 
for  j=l:size(sig,2) 
if  mxm(l,i)==sig(i,j) 
plaoe(i,l)=j; 
end 
end 
end 


ss=0; 

for  i=l;size(sig,l) 
for  j=l:size(sig,2) 
if  sig(i,j)>A 

ifj>=5  &  j<=size(sig,2)-4 
xmtifft(i,[j-4:j+4])=xmtifft(i,[j-4;j+4]).*wo; 
elseif  j<5 

xmtifft(i,[j;j+4])=xmtifft(i,[j;j+4]).*(l-ko*hanning(5)'); 
elseif  j>size(sig,2)-4 
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xmtifft(i,[j-4;j])=xmtifft(i,[j-4:j]).*(l-kc*hanning(5)'); 

end 

end 

end 

end 


BB=1; 

ifBB==l 

B=2;  %BOTTOM  LEVEL 
ss=0; 

for  i=l;size(sig,l) 
for  j=l:size(sig,2) 
if  sig(i,j)<B 

ifj>=5  &  j<=size(sig,2)-4 
xmtifft(i,[j-4;j+4])=xmtifft(i,[j-4;j+4]).*wo2; 
elseif  j<5 

xmtifft(i,[j;j+4])=xmtifft(i,[j;j+4]).*(l+ka*hanning(5)'); 
elseif  j  >size(sig,2)-2 

xmtifft(i,[j-4;j])=xmtifft(i,[j-4:j]).*(l+ka*hanning(5)'); 

end 

end 

end 

end 

end 

%  PEAK  CLIPPING 
elseif  window==2 
CL=12;  %CLIPPING  LEVEL 
xmtifft=xmtifft  *  5  0 ; 
xmtifft=xmtifft( : ,  [  1 7 : 8  0] ) ; 
sig=abs(xmtifft); 
figure(30);plot(sig); 
for  i=l:size(sig,l) 
for  i=l:size(sig,2) 
if  sig(i,j)>CL 
xmtifft(i,j)=CL; 
end 
end 
end 

end 

if  window~=0 
%  insertion  guard  interval 
Ng=16; 

xmtifft=  [xmtifft( :  ,N -Ng+ 1  :N)  xmtifft] ; 
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%xmtifft_afterGI=size(xmtifft) 

ppp=xmtifft*.02; 

elseif  window==0 

ppp=xmtifft; 

end 

%  PAPR  CALCULATION 
p  ow=xmtifft .  *  c  onj  (xmtifft) ; 
powl=pow';  %transpose 

maxi=max(powl);  %lind  peak  ins.  power  of  OFDM  signal 
ave=mean(powl);  %find  average  ins.  power  of  OFDM  signal 
pap=maxi./ave;  %PAPR  array  (PAPR  value  of  each  OFDM  symbol) 
max_PAPR=10*logl0(max(pap))  %  maximum  PAPR  in  dB 
ave_PAPR=10*logl0(sum(pap)/size(pap,2))  %average  PAPR  in  dB 
% - 
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RAY  DOP 


%RAY_DOP 

%  TITLE  :  RAYLEIGH  DOPPLER 

% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 

%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  o=ray_dop(s,M,N,es) 

m=0:M-l; 

randn('seed',s+10); 

prl=randn(l,20); 

randn('seed',s+20); 

pim=i*randn(  1 ,20); 

p=prl+pim; 

p=p/(40^.5); 

rand('seed',s+30); 

e=rand(l,20); 

e=es*oos(2*pi*(e-.5)); 

E=exp(i*2*pi*e'*m/N); 

o=p*E; 

% - 
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RMDECODER 


%RMDECODER 
%  TITEE  :  RM  DECODER 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  TEZEREN,Serdar  Umit,  ETJG,  TURKISH  NAVY 

% - 

funetion  [msg_deooded_line]=rmdeooder(Me,r,m); 
k=l;%  #  of  the  rows=k;#  of  the  columns=2^m 
for  i=l:r; 

k=k+prod(l  :m)/prod(l  ;m-i)/prod(l  :i); 
end 

[GM,k]=rmgenmat(r,m); 

[WWl  ]=WW_mat2(GM,r,m,k); 
rrr=rem(size(Me,2),2^m); 
if  rrr~=0; 

%display('ERROR! ! !  The  #  of  the  bits  at  the  eneoder  output  is  not  equal  to  the  #  of  the 
bits  at  the  deeoder  input'); 

Me=[Me  zeros(l,2^m-rrr)]; 
end 

Me=reshape(Me,2^m,length(Me)/2^m)’;%RESHAPE  Me  AS  2^m  COEUMNS 
cv_row=l; 
for  i=2:k; 
oounter=0; 

for  j= 1 :  size(  WW  1,1); 

dpt=dot_pro(GM(i,:),WWl(j,:))i 

if  dpt==l; 

dpt2_oounter=0;flag=0; 
for  z=2:k; 

dpt2=dot_pro(GM(z,:),WWl(j,0)i 

if  dpt2==0; 

dpt2_oounter=dpt2_oounter+ 1 ; 
end 
end 

if  dpt2_eounter==k-2; 

CV(ev_row,:)=WWl(j,:)i 

eheek(  1  ,ev_row)=i; 
eheek(2,cv_row)=j ; 
ev_row=ev_row+ 1 ; 
eounter=eounter+ 1 ; 

GG(l,i)=eounter; 

end 

end 
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end 


dpO=0;dpl=0;row_dd=l;My=zeros(size(Me,l),length(GM)); 

ss=0; 

for  f=l:size(Me,l); 
for  i=2:k; 

for  j= 1  ;GG(  1  ,i);%(2^(m-r)); 

dp=dot_pro(CV(row_dd,:),Me(f,:)); 

DP(i,j)=dp; 
if  dp==l 
dpl=dpl+l; 
elseif  dp==0; 

dpO=dpO+l; 

end 

row_dd=row_dd+ 1 ; 
end 

if  dpl>dpO 
dotpro=l; 
else  %if  dpO>dpl 
dotpro=0; 
end 

dp0=0;dpl=0; 

WWW(l,i-l)=dotpro; 

My(f,:)=xor(My(f,:),(dotpro*GM(i,:))); 

end 

MM(f,:)=xor(My(f,:),Me(f,:)); 
ttl=0;tt0=0; 
for  i=l  ;size(MM,2); 
ifMM(f,i)==0; 

ttO=ttO+l; 
elseif  MM(f,i)==  1 ; 

ttl=ttl+l; 

end 

end 

if  ttl>ttO 
tt=l; 

elseif  ttO>ttl; 
tt=0; 

elseif  ttO==ttl 
tt=l; 
end 

msg_deeoded(f,  1  )=tt; 
for  i=2:k; 

msg_deooded(f,i)=WWW(l,i-l); 

end 


121 


row_dd=l; 

end 

Mc=xor(MM,Me); 

msg_decoded_line=reshape(msg_decoded',l,size(msg_decoded,l)*size(msg_decoded,2)) 
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RMENCODER 


%RMENCODER 
%  TITLE  :  RM  ENCODER 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  TEZEREN,Serdar  Umit,  LTJG,  TURKISH  NAVY 

% - 

funetion  [enod_msg]=rmenooder(r,m,msgl); 

[GM,k]=rmgenmat(r,m); 
len_msg=size(msg  1 ,2); 

rrr=rem(len_msg,k);  %oheok  whether  padding  is  neeessary 
if  rrr>0; 

msg=[msgl  zeros(l,k-rrr)]; 
elseif  rrr==0; 

msg=msgl; 

end 

len=size(msg,2)/k; 

msg=  reshape(msg,k,len)';%make  the  message  in  a  matrix  form 
Mo=msg*GM; 
for  i=l;len; 
for  j=l:2^m; 

if  rem(Mo(i,j),2)==0; 

Mo(i,j)=0; 

elseif  rem(Me(i,j),2)>0; 

Me(i,j)=l; 

end 

end 

end 

Me;  %ENCODED  MESSAGE  MATRIX 

ened_msg=reshape(Me',l,(size(Me,l)*size(Me,2)));%ENCODED  MESSAGE 
% - 


123 


RMGENMAT 


%RMGENMAT 

%  TITLE  ;  RM  ENCODING  MATRIX  GENERATOR 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 

%  Author  ;  TEZEREN,Serdar  Umit,  LTJG,  TURKISH  NAVY 

% - 

funetion  [GM,k]=rmgenmat(r,m); 

k=l;%  #  of  the  rows=k;#  of  the  eolumns=2^m 

for  i=l:r; 

k=k+prod(l  :m)/prod(l  ;m-i)/prod(l  :i); 
end 

[R]=gm_partl(m);  %ereating  the  monomials  of  the  generator  matrix  (l,xl,x2,....,xm) 

row=m+2;  %oreating  the  seeond  part  of  the  generator  matrix  (xlx2,xlx3, . ,xm- 

Ixm) 

for  i=2;m; 

for  j=i+l:m+l; 

R(row,:)=and(R(i,:),RO,:)); 

row=row+l; 

end 

end 

if  k>row-l;  %  cheek  whether  the  third  part  is  necessary 
%  if  necessary, 

for  i=2:m+l;  %creating  the  third  part  of  the  generator  matrix  (xlx2x3,x  1x2x4, ,xm- 

2xm-lxm) 

for  j=i+l:m; 
for  z=j+l;m+l; 

R(row,:)=R(i,:)&RO,:)&R(z,:); 

row=row+l; 

end 

end 

end 

end 

if  k>row-l;%  check  whether  the  fourth  part  is  necessary 
%  if  necessary, 

%creating  the  fourth  part  of  the  generator  matrix  (xlx2x3x4,xlx2x3x5, . ,xm-3xm- 

2xm-lxm) 
for  i=2:m-2; 
for  j=i+l:m-l; 
for  z=j+l;m; 
for  y=z+l:m+l; 

R(row, :  )=R(i,  :)&RG ,  :)&R(z,  :)&R(y , :); 
row=row+l; 
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end 

end 

end 

end 

end 

if  k>row-l;%  eheek  whether  the  fifth  part  is  neeessary 
%  if  neeessary,  add  the  fifth  part 
for  i=2:m-2; 
for  j=i+l;m-l; 
for  z=}+\:m; 
for  y=z+l:m+l; 

R(row,:)=R(i,:)&RO,:)&R(z,:)&R(y,:); 

row=row+l; 

end 

end 

end 

end 

end 

if  k>row-l;%  eheek  whether  the  sixth  part  is  neeessary 
%  if  neeessary,  add  the  sixth  part 
for  i=2:m-3; 
for  j=i+l:m-2; 
for  z=j+l:m-l; 
for  y=z+l;m; 
for  h=m+l; 

R(row,:)=R(i,:)&RG,:)&R(z,:)&R(y,:)&R(h,:); 

row=row+l; 

end 

end 

end 

end 

end 

end 

if  k>row-l 

display('ERROR! ! !  the  generator  matrix  need  to  be  modified  seventh  part  is  neeessary') 
end 

for  i=l;k; 

GM(i,:)=R(i,:); 

end 

GM;%  RM  ENCODING  MATRIX 

% - 
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ROTM 


%ROTM 

%  TITLE  ;  ROTATE  VECTOR 

% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 

%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  [vp,vn]=rotm(v,m) 

E=length(v); 

m=rem(m,E); 

ii=(l:E)-l; 

isp=rem(ii-m+E,E)+ 1 ; 
isn=rem(ii+m+E,E)+ 1 ; 
vp=v(isp); 
vn=v(isn); 

% - 
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TDA 


%TDA 

%  TITLE  :  EFT  AND  GI  INSERTION  OPERATOR 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Prof  Paul  H.  Moose,  Naval  Postgraduate  Sehool 

%  Modified  by  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

funetion  x=tda(Ng,X) 

[m  N]=size(X); 

%  Perform  inverse  EFT  on  frequeney  values  in  array ,X 
x=ifft(X.’); 

%  Add  preeursor  of  Ng  samples  to  the  beginning  of  the  time  domain  array  for  ehannel 

%  eompensation. 

x=x.'; 

ifNg==0 

x=x; 

else 

x=[x(:,N-Ng+l:N)  x]; 
end 

% - 
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VITERBI 


%VITERBI 

%  TITLE  :  VITERBI  DECODER 


% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  Tan  Kok  Chye,  Naval  Postgraduate  Sehool 

% - 

fune- 

tion[deooder_output,survivor_state,oumulated_metrio]=viterbi(v_G,v_k,ohannel_output) 

v_n=size(v_G,l); 

%  eheek  the  sizes 
if  rem(size(v_G,2),v_k)~=0 
error('Size  of  v_G  and  v_k  do  not  agree') 
end 

if  rem(size(ehannel_output,2),v_n)~=0 
error('ehannel  output  not  of  the  right  size') 
end 

v_L=size(v_G,2)/v_k; 
number_of_states=2^((v_L- 1 )  *  v_k); 

%generate  state  transition  matrix,  output  matrix,  and  input  matrix 
for  v_j=0;number_of_states-l 
for  v_l=0:2^v_k-l 

[next_state,memory_oontents]=nxt_stat(vj,v_l,v_L,v_k); 

input(v j + 1  ,next_state+ 1  )=v_l; 

branoh_output=rem(memory_oontents  *  v_G'  ,2) ; 

nextstate(v  j + 1 ,  v_l+ 1  )=next_state; 

output(vj+l,v_l+l)=bin2deoi(branoh_output); 

end 

end 

state_metrie=zeros(number_of_states,2); 

depth_of_trellis=length(ehannel_output)/v_n; 

ehannel_output_matrix=reshape(ehannel_output,v_n,depth_of_trellis); 
survivor_state=zeros(number_of_states,depth_of_trellis+ 1 ); 

%start  deeoding  of  non-tail  ehannel  outputs 
for  v_i=l  ;depth_of_trellis-v_L+l 
flag=zeros(  1  ,number_of_states); 
if  v_i<=v_L 

step=2^((v_L-v_i)*v_k); 

else 

step=l; 

end 

for  v_j=0:step:number_of_states-l 
forv  l=0:2^v  k-1 
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branch_metric=0; 

binary_output=deci2bin(output(v j+ 1  ,v_l+ 1  ),v_n); 
for  v_ll=l:v_n 

branch_metric=branch_metric+metric(channel_output_matrix(v_ll,v_i),binary_output(v_ 

11)); 

end 

if((state_metric(nextstate(v_j + 1  ,v_l+ 1 )+ 1 ,2)>state_metric(v_j+ 1 , 1  )+branch_metric)|  flag( 
nextstate(v j + 1 ,  v_l+ 1 )+ 1  )==0) 

state_metric(nextstate(v_j + 1  ,v_l+ 1 )+ 1 ,2)=state_metric(v_j +1,1  )+branch_metric; 

survivor_state(nextstate(v  j + 1  ,v_l+ 1 )+ 1 ,  v_i+ 1  )=v J ; 

flag(nextstate(v j + 1 ,  v_l+ 1 )+ 1  )=  1 ; 

end 

end 

end 

state_metric=state_metric(:,2:-l ;  1); 
end 

%start  decoding  of  the  tail  channel  outputs 
for  v_i=depth_of_trellis-v_L+2;depth_of_trellis 
flag=zeros(  1  ,number_of_states); 

last_stop=number_of_states/(2^((v_i-depth_of_trellis+v_L-2)*v_k)); 

for  vj=0:last_stop-l 

branch_metric=0; 

binary_output=deci2bin(output(vj+l,l),v_n); 
for  v_ll=l:v_n 

branch_metric=branch_metric+metric(channel_output_matrix(v_ll,v_i),binary_output(v_ 

11)); 

end 

if((state_metric(nextstate(v_j + 1 , 1 )+ 1 ,2)>state_metric(v_j +1,1  )+branch_metric)  |  flag(next 
state(  V  J +1,1)+1)==0) 

state_metric(nextstate(v_j + 1 , 1 )+ 1 ,2)=state_metric(v_j +1,1  )+branch_metric; 

survivor_state(nextstate(v J + 1 , 1 )+ 1 ,  v_i+ 1  )=v j ; 

flag(nextstate(  V J +1,1)+1)=1; 

end 

end 

state_metric=state_metric(:,2:-l ;  1); 
end 

%generate  the  decoder  output  from  the  optimal  path 
state_sequence=zeros(  1  ,depth_of_trellis+ 1 ); 

state_sequence(  1  ,depth_of_trellis)=survivor_state(  1  ,depth_of_trelhs+ 1 ); 
for  v_i=l  ;depth_of_trellis 
state_sequence(  1  ,depth_of_trellis- 

v_i+ 1  )=survivor_state((state_sequence(  1  ,depth_of_trellis+2-v_i)+ 1  ),depth_of_trellis- 

v_i+2); 

end 

decoder_output_matrix=zeros(v_k,depth_of_trellis-v_L+ 1 ); 
for  v_i=l  ;depth_of_trellis-v_L+l 
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dec_output_deci=input(state_sequence(  1  ,v_i)+ 1  ,state_sequence(  1  ,v_i+ 1 )+ 1 ); 
dec_output_bin=deci2bin(dec_output_deci,v_k); 
decoder_output_matrix(:,v_i)=dec_output_bin(v_k:-l ;  1)'; 
end 

decoder_output=reshape(decoder_output_matrix,l,v_k*(depth_of_trellis-v_L+l)); 
cumulated_metric=state_metric(  1,1); 

% - 
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WWMAT  2 


%WWMAT_2 

%  TITLE  :  CHARACTERISTIC  VECTOR  EINDER 

% - 

%  Thesis  Advisor  :  Prof  M.  Tummala,  Naval  Postgraduate  Sehool 
%  Author  ;  TEZEREN,Serdar  Umit,  ETJG,  TURKISH  NAVY 

% - 

funetion  [WWl]=WW_mat2(GM,r,m,k); 

[GM,k]=rmgenmat(r,m); 
row_d=l;%rows  themselves  and  inverses 
for  i=2;m+l; 

WWl(row_d,:)=GM(i,:); 

WWl(row_d,:)=inverse(GM(i,:)); 
r  0  w_d=r  0  w_d+2 ; 
end 

%double  eombinations 
for  i=2:k; 
for  j=i+l:k; 

WWl(row_d,:)=GM(i,:)  &  GMO,:); 

WWl(row_d+l,:)=inverse(GM(i,:))  &  GM(j,:); 

WWl(row_d+2,:)=GM(i,:)  &  inverse(GM(j,:)); 
WWl(row_d+3,:)=inverse(GM(i,:))  &  inverse(GM(j,:)); 
row_d=row_d+4 ; 
end 
end 

for  00=3  :m-l; 

for  j= 1 :  size(  WW  1,1); 
for  i=oo+l;m+l; 

WWl(row_d,:)=GM(i,;)  &  WWIG,:); 

WWl(row_d+l,:)=inverse(GM(i,:))  &  WWlG,0i 
row_d=row_d+2 ; 
end 
end 
end 

WWl(row_d,:)=ones(l,size(WWl,2));%adding  all  ones  matrix  at  the  end  of  the  matrix 
% - 
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